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S_i ABSTRACT 

Over the last forty years, the 12 C + 12 C fusion reaction has been the subject of 
03 considerable experimental efforts to constrain uncertainties at temperatures relevant 

for stellar nucleosynthesis. Recent studies have indicated that the reaction rate may 
be higher than that currently used in stellar models. In order to investigate the effect 
of an enhanced carbon burning rate on massive star structure and nucleosynthesis, 
I/-) new stellar evolution models and their yields are presented exploring the impact of 

{SJ three different 12 C + 12 C reaction rates. Non-rotating stellar models considering five 

CN different initial masses, 15, 20, 25, 32 and 60Mq, at solar metallicity, were generated 

t-H using the Geneva Stellar Evolution Code (GENEC) and were later post-processed with 

* the NuGrid Multi-zone Post-Processing Network tool (MPPNP). A dynamic nuclear 

t^, reaction network of ~ 1100 isotopes was used to track the s-process nucleosynthesis. 

(S) An enhanced 12 C + 12 C reaction rate causes core carbon burning to be ignited more 

i— I promptly and at lower temperature. This reduces the neutrino losses, which increases 

the core carbon burning lifetime. An increased carbon burning rate also increases 
. — the upper initial mass limit for which a star exhibits a convective carbon core (rather 

than a radiative one). Carbon shell burning is also affected, with fewer convective-shell 
5— i episodes and convection zones that tend to be larger in mass. Consequently, the chance 

of an overlap between the ashes of carbon core burning and the following carbon shell 
convection zones is increased, which can cause a portion of the ashes of carbon core 
burning to be included in the carbon shell. Therefore, during the supernova explosion, 
the ejecta will be enriched by s-process nuclides synthesized from the carbon core s 
process. The yields were used to estimate the weak s-process component in order to 
compare with the solar system abundance distribution. The enhanced rate models 
were found to produce a significant proportion of Kr, Sr, Y, Zr, Mo, Ru, Pd and Cd 
in the weak component, which is primarily the signature of the carbon-core s process. 
Consequently, it is shown that the production of isotopes in the Kr-Sr region can be 
used to constrain the 12 C + 12 C rate using the current branching ratio for a- and 
p-exit channels. 

Key words: nuclear reactions, nucleosynthesis, abundances - stars: abundances - 
stars: evolution 
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1 INTRODUCTION 

Despite the limitations of ID stellar models, their capabil- 
ity to reproduce several observables makes them a funda- 
mental tool to understand stellar nucleosynthesis sites in 
the galaxy. Calculated stellar abundances can be compared 
with observed abundances from meteoritic data or stellar 
spectra. In massive stars (M > 8Mq) the presence of ad- 
vanced burning stages during their evolution and their final 
fate as a supernova explosion provides a useful test-bed for 
many sensitivity studies, which are important to constrain 
uncertainties in input physics. In particular, nuclear reac- 
tion rates are often found to be sources of uncertainty as the 
task of experimentally determining precise cross sections at 
astrophysically relevant energies is often difficult. The 12 C 
+ 12 C reaction is a good example where, despite over four 
decades of research, the reaction rate still carries substantial 
uncertainties because of the nuclear structure and reaction 
dynamics governing the low energy cross section of the fusion 
process ( Strieder|2010 I. The extrapolation of the laboratory 
data into the stellar energy range - Gamow peak energies 
(E ~ 1.5 MeV, or T ~ 0.5 GK) - depends critically on 
a reliable theoretical treatment of the reaction mechanism. 
Present model extrapolations differ by orders of magnitude; 
this affects directly the reaction rate with significant impact 



ditions (other than central density) are less sensitive to an 



on a number of stellar burning scenarios (Gasques et al. 
|2007| ). 

The 12 C + 12 C reaction cross section is characterized 
complex resonance structure associated either with scatter- 
ing states in the nucleon-nucleon potential or wi th quasi- 
molecular states of the compound nucleus 24 Mg (Imanishi 
|1968[ ), which at low energies can be described by a resonant- 
part superimposed on a non-resonant part, where the latter 
is also rather uncertain ( Yakovlev et al.|2010 1. A theory that 
predicts the location and strength of the resonant-part has 
not yet been proposed ( Strieder||2008 1, but resonance char- 
acteristics can be determined either by coupled-channel cal- 
culations or optical model potentials based on, for example, 
a-particle condensates or cluster structures ( Xu et al.|2010 



Betts fc Wuosmaa|1997 and references therein) . Resonances 



have consequently been predicted by both approaches at en- 
ergies ~ 2 MeV (|Michaud fc Vogt||1972| |Perez-Torres et al 



20061 and it was shown that the experimentally observed 



data could be reasonably well reproduced in the framework 
of these models ( Kondo, Matsuse, fc Abe~||1978 1. Yet, none 
of these models provides the quantitative accuracy in res- 
onance parameter predictions, required for a reliable ex- 
trapolation of the data into the stellar energy range. Com- 
plementary to the classical potential model approach, dy- 
namic reaction theories are being developed. They have been 
tested successfully for fusion of spherical nuclei like 16 O + 

R J 1 L 

O (Diaz-Torres et al.|2007[), but the theoretical treatment 



of fusion reactions of two deformed C nuclei requires a 
non-axial symmetric formalism for a fully reliable treatment 
( Diaz-Torres|2008 |. 

Taking a phenomenological approach a resonance with 
strength (loj) ~ 3.4 x 10 -7 eV has been invoked to correct 



the ignition depth of neutron star superbursts ( Cooper et al 



2009 1, which are believed to be caused by ignition of carbon- 
burning reactions, triggering a thermonuclear runaway in 
the crust of a neutron star. Type la supernovae should also 
exhibit changes to the ignition characteristics, but these con- 



enhancement in the carbon burning rate ( Cooper et al. 2009 



Iapichino & Lesaffre 20101. The possible existence of such 
a resonance, associated with a pronounced 12 C + 12 C clus- 
ter structure of the compound nucleus 24 Mg, represents a 
source of uncertainty. 

Alternatively, the reaction rate may not be dominated 
by resonances at lower energies because of predictions that 
the cross section drops much steeper than usually antici- 
pated due to a fusion hindrance reported in heavy-ion reac- 



tions (see for example, |Jiang et al. |2004 2007). The con- 
sequences of the hindrance phenomenon for the 12 C+ 12 C 
reaction in astrophysical scenarios was examined by |Gasques| 
(20071, where it was demonstrated that hindrance is 



et al. 



much more significant in the pycnonuclear regime than the 
thermonuclear regime, but does exhibit a noticeable effect 
on the yields of massive stars. The reduced rate, by approx- 
imately a factor of 10-100 at carbon burning temperatures 
(see their Fig. 1), increases the temperature with which car- 
bon burning occurs and therefore affects the nucleosynthesis. 
Changes in the yields were generally rather small, but some 
specific isotopes, such as 26 Al, 40 Ca, 46 Ca, 46 Ti, 50 Cr, 60 Fe, 
74 Se, 78 Kr and 84 Sr, exhibited larger changes most likely due 
to the increased neutron density exhibited by the burning 
of neutron sources at higher temperatures. 

The wide range of presently discussed model predictions 
requires new experimental effort to reduce the uncertainty 
range. However, the measurements towards low energies are 
extremely difficult, because the low cross section (a <C 1 
nbarn) limits the experimental yield to an event rate be- 
low the natural and beam induced background events in 
the detectors. Particle measurements are difficult because of 
the limited energy resolution of the particle detectors which 
makes a separation of the particle groups extremely difficult 
at the low count rate conditions. Beam induced background 
from reactions on target impurities is therefore difficult to 



distinguish from the actual reaction products (Zickcfoosc 
II2010D . The measurement of secondary gamma radia- 



let al. 



tion associated with the particle decay is also handicapped 
by natural and cosmic ray induced background radiation 
(Strieder 20101. While recent experiments suggest an in- 



crease in the low energy S-factor indicating the possibility 



of narrow resonances at lower energies ( Barron-Palos et al 



2006) |Aguilera et al.|[2006] |Spillane et ai.||2007[ ), the confir 



mation of the results and the experimental pursuit towards 
lower energies is stalled due to the present inability to dif- 
ferentiate the reaction data from the different background 
components JZickefoose et al. 20101. Improved experimen- 



tal conditions requires the preparation of ultra-pure target 
materials for experiments in an cosmic ray shielded under- 
ground environment ( Strieder|2010 1. 

The three dominant carbon burning reactions, with Q- 
values, are 



12 C( 12 C,a) 20 Ne, 
12 C( 12 C,p) 23 Na, 
12 C( 12 C,n) 23 Mg, 



= +4.617 (1.1) 
= +2.240 (1.2) 
= -2.599. (1.3) 



During carbon-burning, the a- and p-channels dominate 
with the n-chan nel making up less t han 1 per cent of all 12 C 
+ 12 C reactions ( Dayras et al.|l977 l. At this stage, the com- 



position of the star is largely C and O, with the initial 
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ratio of 12 C to 16 at this stage largely governed by the 
12 C (a, 7) 16 O reactions occurring during helium-core burn- 
ing. Carbon-core burning occurs at a central temperature 
~ 0.7 GK and produces mainly 20 Ne and 24 Mg, since ~ 99 
per cent of 23 Na synthesised through the p-channel is de- 
stroyed via efficient 23 Na(p, q) 20 Ne a nd 23 Na(p, 7) 24 Mg re- 
actions ( jArnett fc Thielemann|1985[ ). Carbon-core burning, 



suiting in an s-process with a higher neutron density and a 



which is convective for stars with initial mass M < 20A/ ( 
and radiative for M > 20A1q (see for example 



G 



Hirschi 



et al.|[2005 1, is followed by convective carbon-shell burning 
episodes at temperatures ~ 0.8 — 1.4 GK. The number of 
episodes and the spatial extent of each shell differs between 
massive stars of different initial mass as the development of 
the carbon shells is sensitive to the spatial 12 C profile at 
the end of helium-core burning; the formation of a convec- 
tive carbon-shell often lies at the same spatial coordinate 
as the top of the previous convective shell ( Arnett|1972 El 
Eid et al.p 004). The presence of a convective carbon core 
depends on the CO core mass as both the neutrino losses 
and energy generation rate depend on the density, which de- 
creases with increasing CO core mass ( |Arnett|1972l|Woosley 
& Weaver 1986; Limongi et al. 20001. Consequently, mech- 
anisms that affect the CO core mass or the carbon burning 
energy budget, such as rotation ( jHirschi et al.||2004[) and 
the 12 C abundance following helium burning (Imbriani et al. 
2001||E1 Eid et al. 12009] ), will affect the limiting mass for the 
presence of a convective core. 

Massive stars are a site for the s process, which starts 
during helium-core burning and also occurs during the fol- 
lowing carbon burning stages. S-process nucleosynthesis also 
occurs in the helium-shell via the 22 Ne neutron source, but 
this process is marginal compared to the s process operat- 
ing in the helium-core or the carbon shells (see for example 
The et aL]|2007 |. Beyond carbon burning, the temperature 
becomes high enough in the interior (~ 2 GK) for photodis- 
integration reactions to destroy heavy nuclides. Because the 
s process can probably occur during both central and shell 
carbon-burning, one can expect that changes in the 12 C + 
12 C rate affect the stellar structure and nucleosynthesis and 
therefore also the s process. 

The 22 Ne neutron source, which is formed during he- 
lium burning via the 14 N(q, 7) 18 F(/3+) 18 Q(q, 7) 22 Ne reac- 
tion chain is the main neutron source flPeters||l968l |Couch| 



et al. 1974 Lamb et al. 19771. As the temperature ap- 
proaches 0.25 GK near the end of helium-burning, 22 Ne(a, 
n) 25 Mg reactions become efficient (Busso & Gallino 



1985 



Raiteri et al. 19911. During this phase a 25Mq star, for 



example, has a neutron density n n ~ 10 B cm" 3 a nd a neu- 
tron exposu re r n ~ 0.2 mb" 1 (see for instance Pignatari 

The 22 Ne source be- 



et al. 



2010 



and references therein) 
comes efficient in a convective environment and heavy ele- 
ments formed through neutron captures are mixed out from 
the centre of the star. Some of these abundances will be mod- 
ified by further explosive nucleosynthesis later in the evolu- 
tion, but will otherwise survive long enough to be present 
in the supernova ejecta and contribute to the total yields 
of the star. Consequently, 22 Ne in massive stars is the dom- 
inant neutron source responsible for the classical weak-s- 



process component (iTruran & Iben 1977 Prantzos et al 



1987| |Kappeler et al.|1989||Raiteri et al.|1991 1 

Any remaining Ne present at the end of helium-core 
burning is later reignited during carbon-shell burning re- 



lower neutron exposure ( 



10 1 



cm and r„ 



0.06 



mb~ ; Raiteri et al. 1991 1. The increased neutron density 
is responsible for changing the branching ratios of unsta- 
ble isotopes, which is particularly important for branching 
isotopes, such as 69 Zn, 79 Se and 85 Kr, since they inhabit 
positions in the isotope chart of nuclides where different 
s-process paths across the valley of stability are available 
( Kappeler et al.|1989| ). The increase in neutron density is re- 
sponsible for opening the s-process path so that the carbon- 
shell burning contribution to specific isotopes, such as 70 Zn, 
86 Kr and 80 Se, may be relevant (see for example Raiteri et al. 
IMll|The~et al.|2007 l. 

Another potential neutron source is 13 C, whic h is 
formed through the 12 C(p,7) 13 N(/3 + ) 13 C reaction chain |Ar~ 



nett fc Truran|1969 1. During carbon-core burning this neu- 
tron source, via the 13 C(a,n) 16 reaction, becomes efficient 
which results in an s-process in the c arbon-core with a typ 



ical neutron density n„ = W 7 cm 3 (Arnett & Thielemann 



1998 1. The abundance of 13 C is dependent 



[19851 iChieffiet al. 
on the N(7,p) C reaction, which dominates the depletion 
of 13 N at temperatures above 0.8 GK. The 22 Ne neutron 
source is the dominant neutron source when the tempera- 
ture rises above such a temperature, although the 13 C neu- 
tron source may also provide an important contribution to 



the total neutron exposure (Clayton 1968 Arcoragi et al 



1991 1. In any case, the carbon-core s process occurs primar- 



ily in radiative conditions with a relatively small neutron ex- 
posure and any heavy elements synthesised via the ensuing 
neutron-captures usually remain in the core (see however the 
discussion on overlapping convection zones in sQ; photodis- 
integration and the supernova explosion process will ensure 
that these elements are not present in the final ejecta and do 
not contribute to the final yields of the star (see for example, 
Chieffi et al.|l998 |. 



A preliminary study ( Bennett et al.pOlO i) found that 
changes to the total 12 C + C rates within a factor of 10 
affect the convection zone structure and nucleosynthesis of a 
25 Mq star at solar metallicity. The main conclusions were 
an increase in the carbon-burning shell contribution to the 
s-process abundances by two different scenarios. The first, 
applicable to the case where the rate was increased by a fac- 
tor of 10, was due to the presence of large carbon-burning 
shells that 'overlapped'. In this situation, the second carbon- 
burning shell was polluted with ashes from the first carbon- 
burning shell, modifying the overall composition. The second 
scenario, applicable to the case where the rate was reduced 
by a factor of 10, was an increase in neutron density asso- 
ciated with the neutron source, 22 Ne, burning at a higher 
temperature in the convective shell. The overall increase in 
the abundances of most isotopes with 60 < A < 90 was 
approximately 0.1 to 0.4 dex. Strongly enhanced rates were 



also investigated (Bennett et al. 2010 3), which show that 
the presence of a larger convective core has a significant im- 
pact on the total yields, since the convective core adds an 
additional neutron exposure towards the total contribution 
of s-process yields; abundances of many heavy nuclides in- 
creased by up to ~ 2 dex. However, no comparison could be 
made with observations as a 25 Mq stellar model (at solar 
metallicity) was the only one considered. 

In this paper, a sensitivity study is made over a set 
of massive star models, at solar metallicity, to determine 
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whether a comparison between the yields and the solar sys- 
tem abundances can constrain the 12 C + 12 C rate, ^ex- 
plains the models and the choice of input physics in the sim- 
ulations. In Sj3] the changes in stellar structure are analysed. 
Sj4] describes the nucleosynthesis, focusing on the s process 
during carbon-core and carbon-shell burning. ij5] presents 
the yields. The discussion and conclusions can be found in 
!|6]and ^respectively. 



2 COMPUTATIONAL APPROACH 



2.1 The 12 C + 12 C reaction rates 



We build on the previous work ( Bennett et al.|2010 j) where 
three carbon burning rates in a 25M q star were considered. 
These are the Caughlan & Fowler (|1988[) 'standard' rate 



(ST) and two enhanced rates: an 'upper limit' rate (CU) and 
an intermediate rate (CI), the latter of which is a geometric 
mean of the ST and CU rates. The CU rate is the ST rate 
including a reasonance of strength (wy) = 6.8 x 10 -5 eV at 
a centre-of-mass energy -Z5 C om = 1.5 MeV. This choice of res- 
onance originates from a preliminary particle spectroscopy 
experiment on 12 C + 12 C obtained at the CIRCE radioac- 
tive beam facility in Caserta/Napoli, Italy ( |Terrasi et al.| 
2007). Although the CI rate was determined via a geometric 



mean, a resonance that would replicate the peak at 1.5 MeV 
for this rate would have a magnitude of (ury) — 3.4 x 10~ 7 
eV. The top panel of Fig. [l] shows the Maxwellian-averaged 
cross-sections of the reaction rates as a function of temper- 
ature. The bottom panel shows the reaction rates relative 
to the ST rate. As indicated by Fig. [T] the peak of the CU 
and CI rates is at ~ 0.5 GK and is a factor of approximately 
50, 000 and 250 times the ST rate at that temperature re- 
spectively. The choice of branching ratio for the a- and p- 
exit channels is 13:7, which is valid within the energy range 
4.42 < E com < 6.48 MeV ( jAguilera et~ll]p006| . It is as- 
sumed in this work that the branching ratio is preserved to 
lower centre of mass energies. For the n-exit channel, we use 
the branching ratio from |Dayras et al.| ( |1977[ ) . 



2.2 Stellar models 

Non-rotating stellar models at solar metallicity (Z=0.02) 
were generated using the Geneva Stellar Evolution Code 
(GENEC), with a small nuclear reaction network that takes 
into account the reactions important for energy generation. 
Five masses were considered for each carbon-burning rate, 
which are 15, 20, 25, 32 and 60 Mq, for a total of 15 stellar 
models. These will be referred to as XXYY where XX is the 
initial mass of the star in solar masses and YY denotes the 
rate and is 'ST', 'CI' or 'CU' for the standard, intermedi- 
ate and upper limit rates respectively. The reason for this 
choice of initial masses is to provide yields data over a range 
of masses with approximately even spacing in log-space. 

GENEC is described in detail in Eggcnbcr ger et al.| 
12008), but some important features are recalled here for 
convenience. The Schwarzschild criterion for convection is 



used and convective mixing is treated as a diffusive process 
from oxygen burning onwards. No overshooting is included 
except for hydrogen- and helium-burning cores, where an 
overshooting parameter of a — 0.2Hp is used. Neutrino loss 
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Figure 1. Top panel : Maxwellian-averag ed cross-sections for 12 C 
+ 12 C rates used in ( jBennett et al. |2010| 3) and also in this study. 
The three rates are the Caughlan & Fowler (1988) 'standard' rate 
(ST), an upper limit rate (CU) and an intermediate rate (CI). 
The CI rate is a geometric mean of the ST and CU rates. Bottom 
panel: The Maxwellian-averaged cross-sections relative to the ST 
rate. 



rates are calculated using fitting formulae from Itoh et al. 
(119891), which are the same as those of the more recent evalu- 



ation from Itoh et al. ( 1996 1 for pair and photoneutrino pro- 
cesses. The initial abundances used were those of lGrevesse fcl 
Noels (1993), which correspond directly to the OPAL opac- 
ity tables used (Rog ers et al.|1 996). For lower temperatures, 
opacities from |Ferguson et al.| ( |2005| ) are used. 

Several mass loss rates are used depending on the ef- 
fective temperature, T e g, and the evolutionary stage of the 
star. For main sequence massive stars, where logT^ > 3.9, 



mass loss rates are taken from Vink et al. (2001 1. Otherwise 
the rates are taken from de Jager et al. ( 1988 1. However, for 



lower temperatures (logT e ff < 3.7), a scaling law of the form 



M 



-1.479 x 10~ 



-0 



(2.1) 



is used, where M is the mass loss rate in solar masses per 
year, L is the total luminosity and Lq is the solar luminos- 
ity. For a recent discussion on mass loss rates in the red- 



supergiant phase, see Mauron & Jossclin (2011 1. During the 



Wolf-Rayet (WR) phase, mass loss rates by Nugis & Lamers 
((20001) are used. 

In GENEC the reaction rates are chosen to be those 



of the NACRE compilation; Angulo et al. ( 1999 1 for the 
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experimental rates and from their websiteQ for theoretical 
rates. However, there are a fe w exceptions. The rate of 
Mukhamedzhanov et al] j2003[) was used for 14 N(p, 7) 15 



below 0.1 GK and the lower limit NACRE rate was used 
for temperatures above 0.1 GK. This combined rate is very 



similar to the more recent LUNA rate (Imbriani et al 



2005) at relevant temperatures. The Fynbo 



used for the 3a reaction and the 
was used for 12 C(a,7) 16 Q. The 



Kunz et al 



(2005) rate was 



(2002) rate 



Ne(a,n) 2a Mg rate was 
taken from Jaeger et al. (20011 and used for the avail- 
able temperature range (T ^ 1 GK). Above this range, 
the NACRE rate was used. The 22 Ne(a,n) 25 Mg rate com- 
petes with 22 Ne(cv, 7) 26 Mg for a— particles. For this rate, 
the NACRE rate was used. The 16 O neutron poison is ef- 
fective at capturing neutrons, forming 17 0, which can ei- 



ther resupply the 'recycled' neutrons via the 17 O(a,n) 20 Ne 
reaction or undergo the competing reaction 17 0(a, 7) 21 Ne. 
For 17 0(a, n) 20 Ne the NACRE reaction is use d and for the 
17 0(a, 7) 21 Ne reaction the correction of the Caughlan & 



Fowler (1988) rate by Descouvemont (19931 is applied 



The models were calculated for as far into the evolution 
as possible, which for most models is after or during the 
silicon-burning stage. The models that ceased before silicon 
burning were the 15CI, 15CU, 60CI and 60CU models, which 
proceeded to oxygen-shell burning, and the 20CI and 20CU 
models, which proceeded to just after the oxygen-shell burn- 
ing stage. The s-process yields are not significantly affected 
by hydrostatic burning stages following oxygen burning be- 
cause most of the isotopes produced via the s process will be 
destroyed by photodisintegration and the choice of remnant 
mass for the supernova explosion, which defines the bound- 
ary between matter that falls back onto the remnant and 
matter that forms supernova ejecta, reduces the impact of 
nucleosynthesis that neon, oxygen and silicon burning stages 



would have on the total yields (see also E 5.1 1 . However, it 



must be noted that there will be explosive burning processes 
during the supernova explosion and photodisintegration oc- 
curring at the bottom of the convective carbon, neon and 
oxygen shells during the advanced stages, which will affect 
the abundances (see for example Rauscher et al.|[2002 Tur 
et al.|2 009). In this work the contribution of explosive burn- 
ing and photodisintegration to the total yields is not con- 
sidered. 

Since the 12 C + 12 C reactions do not become efficient 
until after helium-core burning, the CU and CI models for 
a particular choice of initial mass were started just before 
the end of helium-core burning using the ST model data 
as initial conditions, reducing some of the computational 
expense. 



the Newton-Raphson scheme, with the output temperature, 
density and the distribution of convection (and radiation) 
zones from GENEC as input. Additional features have been 
included to enhance the calculations or save on unnecessary 
computations. Sub-timesteps are inserted where appropriate 
to improve convergence in the case where the timescale of 
reactions is smaller than the stellar evolution timestep. Also, 
the nuclear network is dynamic, adding or removing isotopes 
from the network depending on the stellar conditions (up to 
the maximal network defined in Table [TJ . This is useful in 
reducing the number of computations associated with nu- 
clear reactions where the change in abundance is zero or 
negligible. The same (adaptive) mesh used in GENEC was 
used for the post-processing calculations. 

The nuclear networks used are shown in Fig. [2] The iso- 
topes used in each network are discriminated depending on 
whether they are involved in reactions important for energy 
generation (featured in both the stellar model and the post- 
processing tool) or not (featured only in the post-processing 
tool). GENEC uses a skeleton network of 31 isotopes, which 
is the same network used in previous GENEC models (see 
for example |Hirschi et al.||2004] [2005). This network is a 
combination of fundamental isotopes relevant for pp-chain 
reactions, the CNO tricycle and helium burning and a net- 
work similar to the a 7 network of |Hix et al.| ( | 1998p , enacted 
during the advanced burning stages, which reduces the com- 
putational expense associated with a larger network without 
causing significant errors in energy generation rates. The 
isotopes included in the network for MPPNP are specified 
m Tableland are shown in Fig. [2] Five isomeric states are 
also included, which are treated as separate nuclei from their 
ground state equivalents. These are 26 Al m , S5 Kr m , 115 Cd m , 
176 Lu m and 180 Ta m 

The reaction rates in MPPNP were set to those used 
in the skeleton network of GENEC, as specified in j ]2.2[ 
for the same reactions. Additional reactions are taken from 
the default setup of MPPNP and are specified as follows: 



charged particle reactions are from Angulo et al. ( 1999 1 and 



Iliadis et al. (2001 1. /3-decays and electron captures are from 



Oda et al. (19941, Fuller et al. (19851 and Aikawa et al 



(20051. Neutron captures are from the Karlsruhe astrophys- 



ical database of nucleosynthesis in stars (KADoNiS) (Dill 



mann et al. 20061. For reactions not found in these refer- 



ences, reaction rates from the Reaclib database]^] were used, 
which incorporates a compilation of experimental rates and 



theoretical rates from NON-SMOKER (Rauscher & Thiele- 



mann 2000 20011 



2.3 Post-processing 

The NuGricj^] Multi-Zone Post-Processing tool (the parallel 
variant; MPPNP) is described in |Herwig et al. | ( |2008[ ) and 
Pignatari et al. (2011, in prep.). See also appendix [A| for de- 
tails of the parallel implementation. The system of equations 
for the rate of change of abundances of isotopes is solved 
using an implicit finite differencing method combined with 



3 STELLAR STRUCTURE AND EVOLUTION 

3.1 Hydrogen and helium burning 

The evolution of each stellar model during hydrogen- and 
helium-burning is given entirely by the ST models, as the 
CI and CU models were started using the profile just before 
the end of helium burning. Figure [3] shows the Hertzsprung- 
Russell (HR) diagram for all models, which shows that the 
evolutionary tracks for all models follow their course in the 



1 http:/ /pntpm3. ulb.ac.be/Nacre/nacre. htm 

2 http:/ /forum. astro. keele.ac.uk:8080/nugrid 



3 http:/ /nucastro.org/reaclib.html 
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Neutron number 



Figure 2. Chart of isotopes indicating the nuclear reaction networks used in this work: GENEC (blue squares) and MPPNP (pale red 
squares). The network used by MPPNP includes all stable isotopes, which are indicated by black squares. The outer boundary to each 
side of the valley of stability indicates the position of all currently known isotopes, including heavy transuranic isotopes. Parallel grid 
lines indicate values of Z or N that are magic as specified in the nuclear shell model (2,8,20,28,50,82,126). 



HR diagram primarily during the hydrogen- and helium- 
burning phases and are not modified by enhanced rates. The 
reason for this is that the surface evolution of the stellar 
models is unaffected by changes in the carbon-burning rate, 
which is a consequence of the small timescale for burning 
associated with advanced burning stages in massive stars; 
the envelope has insufficient time to react significantly to 
changes in core properties. 

Overall, the ST models are very similar to those pre- 
viously published by the Geneva group, such as the non- 



rotating stars of Meynet & Macdcr ( 2003 1 and Hirschi et al 
( 2004[ ). The 15, 20 and 25 Mq model stars evolve to- 
wards the red and remain as red supergiants (RSGs) dur- 
ing the advanced stages of evolution. The 32 and 60 Mq 
model stars evolve towards the Humphreys-Davidson limit 
at logT c ff ~ 3.8 before becoming WR stars 

The 32 Mq proceeds to the WR phase during helium- 
burning. This is because the mass loss is strong enough for 
the star to expel the entire hydrogen envelope during helium- 
burning, with the composition of the remaining envelope rich 
in helium. The lower opacity of the helium-rich envelope 
lowers the radius and favours evolution towards the blue 
27.3.2). 



(Maeder 2009 



The deviations from the ST track 
for the CI or CU tracks for this mass are slightly larger than 
for other masses. These deviations are generally of the order 
of 0.1 per cent with a maximum deviation of 0.01 in logTcff 
(~ 2 per cent), which occurs during the rapid transit to the 
blue after helium burning. 



phase, the star enters the first 'loop' towards the blue (at 
logT c ff — 4.4), which occurs because of mass loss being high 
enough to expose the helium-rich outer layer. Following the 
first loop to the blue, helium-burning is ignited. During this 
phase the core shrinks, lowering the core fraction, q, favour- 



The 60 M Q 
hydrogen-burning. 



star becomes a WR star just after 
At the end of the hydrogen-burning 



ing evolution to the red ( Maeder||2009| §27.3.2). However, 
the star approaches the Humphreys-Davidson limit in the 
HR diagram during the evolution and the mass loss becomes 
high enough to, eventually, peel away the envelope, expos- 
ing the helium-burning core (q ~ 75 per cent during helium- 
burning). The star consequently evolves towards the blue 
(at logTcff ^ 5.0). 

3.2 Carbon burning 

Unlike the surface evolution, the interior evolution of the 
star is modified significantly by the enhanced carbon burn- 
ing rates and changes to the central evolution of the star are 
important in order to assess changes to the main burning 
regimes. 

Figure [4] shows T c — p c diagrams for the 15, 20 and 25 
Mq models, separated into panels by initial mass. The en- 
hanced rate models in all cases (including the 32 and 60 
Mq models) ignite carbon burning at lower temperatures 
and densities, which consequently affects the evolution of 
the central properties of the star. This is seen, for example, 
in the top and middle panels of Fig. [4] where the curves for 
the CI and CU cases deviate away from that of the ST case 
towards the higher temperature (at a given density) side of 
the curve (see also column 7. in Table |2J. The tendency to 
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Figure 3. Hcrtzsprung-Russell diagram for all models. Solid black lines refer to ST model tracks, dashed blue lines refer to CI model 
tracks and the dotted red lines refer to the CU model tracks. The tracks indicate that the enhanced rates do not affect the surface 
evolution, since changes in the carbon-burning rate do not affect the surface properties. The tracks exhibited by the 32 Mq and 60 Mq 
models show evolution into the WR phase, which is explained by mass loss. 



deviate in this direction is caused by the presence of a con- 
vective core. This is verified in the bottom panel for the case 
of the 25CU model whereby the 'kink' at carbon ignition is 
larger than that of the 25ST and 25CI models, since the CU 
model is the only 25 Mq model to have a convective core 
(see also Fig. [7|. 

Figure [4] shows the impact that the enhanced carbon 
burning rates have on the central evolution during carbon 
burning. However, despite the deviations, many of the mod- 
els at a particular mass are similar, especially the 25 Mq 
models. Figure [5] shows T c — p c diagrams for the 32 and 
60 Mq, which are also quite similar. In the case of Fig. 
[5| the 32 and 60 Mq models exhibit significant mass loss 
during the hydrogen- and helium-burning stages such that 
the total mass during the advanced burning stages is very 
similar (~ 13Mq). Combined with the fact that the helium 
cores at this stage are qualitatively similar, the models from 
this point onwards evolve similarly, with the 32CI and 60CI 
models entering the more degenerate region of the diagram. 
Consequently, the tracks follow similar paths dependent on 
the choice of 12 C + 12 C reaction rate. 

Kippenhahn diagrams for all models are presented in 
Fig. [6j [7| and [8] with the shaded regions corresponding to 
convection zones and the intermediate regions correspond- 
ing to radiative zones. The total mass is given by the thin 
black line at the top of each diagram. Overall, Fig. [6] [7] and 
[8] show that the convection zone structure of the carbon- 
burning stage is heavily modified by the increased rates, 
particularly for the CU cases where a convective carbon- 
core is present over the entire mass range considered. The 



presence of a convective carbon-core is important for nucle- 
osynthesis as the convective mixing provides more fuel for 
carbon-burning and the carbon-core s process. The mass loss 
increases significantly with initial mass, but does not change 
much with the 12 C + 12 C rate. Small deviations in the mass 
loss, which are less than 1 per cent, are due to the increased 
lifetime of the core carbon burning stage in the CI and CU 
models (see Table [If. 

Model data complementary to Figures |6] [7] and H] are 
presented in Table [2] which specify properties pertain- 
ing to convection zones during carbon burning. Column 2. 
('Core/Shell') identifies the presence, or not, of a convec- 
tive core or shell and labels the shells in chronological order 
during the evolution. The other columns specify the lifetime 
of the convection zon^] (tc) in years, the lower and upper 
limits in mass coordinate of the convection zone (Mi ow and 
M upp respectively, in Mq), the size of the convection zone 
in mass (AM, in Mq) and the temperature (T, in GK), 
density (p, in g cm -3 ) and the mass fraction abundances of 
12 C and 16 [Xi2 C and Xibq respectively) at the onset of 
convection at position Mi ow . 



4 Many of the convective shells persist until the presupernova 
stage. In models 15CI, 20CI, 25ST, 25CI, 25CU, 32CI, 32CU and 
60CU however, the carbon shell shrinks because of the influence 
of another burning stage (such as neon or oxygen burning) . The 
convective carbon shell can therefore feature a rather complicated 
structure through the following advanced stages. In these cases, 
the lifetime is calculated from the onset of convection to the point 
where the convective shell shrinks significantly in size. 
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Table 1. Nuclides included in the nuclear reaction network used 
for the post-processing calculations 




Figure 4. T c — p c diagram for all 15 (top panel), 20 (middle panel) 
and 25 Mq (bottom panel) models. The straight line in each panel 
indicates the location in the diagram where the ideal gas pressure 
is equal to the electron degeneracy pressure; P ga s = P e ,deg- Ig- 
nition points for convective core carbon burning are indicated by 
the annotation. 




6 7 
l°g(p c [g ctoT 



Figure 5. T c — p c diagram for all 32 and 60 Mq models. The 
straight line indicates the location in the diagram where the ideal 
gas pressure is equal to the electron degeneracy pressure; Pgas = 

Pe,deg- 
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presence of a convective carbon core with a value between 
20 and 25 Mq, which is consistent with previous models 



(Heger et al. 2000 Hirschi et al. 20041. For model 25CI 



a strong convective shell is ignited slightly off-centre (at a 
mass coordinate of 0.436Mq) and model 25CU exhibits a 
large convective carbon core. In all CU models the carbon- 
core burning stage is convective, which, in models 25CU, 
32CU and 60CU, replaces the radiative cores. In Model 25CI 
the first carbon shell ignites close to the centre and mod- 
els 20CI and 15CI have larger convective cores. Considering 
these facts and the presence of a convective core in every 
CU model, one can hypothesise that the limiting mass for 
the presence of a convective carbon core increases with the 
carbon burning rate, which will consequently represent a 
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logio(Time left until collapse) [yr] logiu(Time left until collapse) [yr] 

Figure 6. Kippenhahn diagrams for the ST, CI and CU models for initial masses of 15 and 20 Mq. Shaded regions correspond to 
convection zones. The major central burning regimes are indicated by the text. 



source of uncertainty for the presence of a convective core 
near to the limiting mass of ~ 22Mq . A firm verification of 
the limiting mass for the CI case would however require a 
finer grid of stellar models between 20 and 25 Mq. 

The sizes, in mass, of the carbon-burning zones (column 
6 in Table [2| are generally larger in CI and CU models. This 
affects the 2 C abundance profile within the star and con- 
sequently the number of carbon-burning shells during the 



evolution. The Kippenhahn diagrams for the 15 and 20 Mq 
models (Fig. [6| demonstrate this effect well; the 15ST and 
20ST models have many carbon burning shells where the 
ignition of a successive shell lies at a position that corre- 
sponds to the maximum coordinate reached by the previous 
convection zone. 

As the rate is increased, the tendency for convective 
shells to 'overlap' (where the lower bound in mass of the 
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Figure 7. Kippenhahn diagrams for the ST, CI and CU models for initial masses of 25 and 32 Mq. Shaded regions correspond to 
convection zones. The major central burning regimes are indicated by the text. 



convective region extends below the upper bound of the pre- 
vious convection zone) is increased. All CU models, except 
the 15CU model, show this overlap, which occurs between a 
convective carbon core and the first convective carbon shell. 
The amount of overlap between the carbon core and the 
first carbon shell, and the first and second carbon shells, in 
the 20CI model (in Fig. [6| is also much larger than that 



in the 20ST model. This overlap effect occurs because suc- 
cessive carbon-shell burning episodes, caused by ignition of 
residual 12 C fuel left over from previous burning stages, can 
occur at a lower temperature and density or with a lower 
abundance of 12 C fuel (see column 9 of Ta ble [2| . This effect 
has been noted previously by Chieffi et al. ( |1998[ ) and in the 
preliminary studies (Be nnett et al.|2010^ ,b). 
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Figure 8. Kippcnhahn diagrams for the ST, CI and CU models 
for initial masses of 15 and 20 Mq. Shaded regions correspond 
to convection zones. The major central burning regimes are indi- 
cated by the text. 



The total energy generation of the 12 C + 12 C reaction 



is given by ( |Woosley et al.|2002[ | 



e„uc( 12 C) « 4.8 x 10 18 F 2 ( 12 C)pAi2,i 2 erg 



s" 1 (3.1) 

where V" 2 ( 12 C), is the number abundance of 12 C (Y — 
X/A), p is the density and Ai2,i2 is the nuclear reaction rate, 
which is dependent on temperature. For a given density and 
abundance, an increased 12 C + 12 C rate increases the en- 
ergy generation rate from nuclear reactions. The effect this 
has on the ignition conditions (temperature and density) for 




Carbon burning rate 

Figure 9. Ignition temperatures for core carbon burning for all 
models. 




Carbon burning rate 

Figure 10. Ignition densities for core carbon burning for all mod- 
els. 



core carbon burning are displayed in Fig. [9] and 10 (the ig- 
nition point is defined as the point in time when the central 
mass fraction abundance of 12 C is 0.3 per cent lower than its 
maximum value). An increased rate allows a star to reach 
the required energy output to support the star against grav- 
itational contraction at a lower temperature (and also lower 
density). Note also the dependence on initial mass, with ig- 
nition conditions favouring higher temperatures and lower 
densities with increasing initial mass. In the case of lower 
ignition temperatures and densities, the convective core ig- 
nites more promptly in the CI and CU models. Changes to 
the ignition conditions and the 12 C abundance at the start 
of core carbon burning are responsible for the increased like- 
lihood of having overlapping convection zones. 

The lifetime of convection zones is generally longer in 
the CI and CU models. This could be perceived as counter- 
intuitive, since with an enhanced rate one would expect that 
the 12 C fuel would be expended more rapidly. However, the 
burning takes place in lower temperature and density con- 
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Table 2. Stellar structure properties for carbon-burning cores and shells at the onset of convection. Shells are labelled in chronological 
order. T conv is the lifetime of the convection zone, Mj ow and M upp are lower and upper mass coordinates for the location of the zone. 
AM is the size of the zone in mass, T and p are the temperature and density of the zone at Mi ow and Xi2q and Xi6q are the 12 C and 
16 O mass fraction abundances within the convection zone. 
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ditions, which affect the neutrino losses. Table [3] shows the 
energy generation terms for nuclear reactions (e nU c) and neu- 
trino losses (e„) at the centre of the star when the mass 
fraction of 12 C is half the amount available just prior to 
carbon-core burning. The proportion of neutrinos formed 
by various neutrino processes are also specified in Table |3j 
which are given as fractions, /, of the total neutrino losses 
(in per cent). These processes are pair production (/ pa ir), 
photoneutrino interactions (/ p hot) and the rest (/ rC st), which 
are bremsstrahlung, recombination and plasmon decay pro- 
cesses ( Itoh et al.||1996| ). Neutrino formation through these 
last three processes is negligibly small at carbon burning 
temperatures. 

As shown by Table [3] the energy generation rate from 
nuclear reactions and the neutrino losses are reduced in the 
CI and CU models, although an increase in energy gener- 
ation rate is seen in models 25CU, 32CU and 60CU from 
their CI counterparts. This increase is due to the presence 
of the convective carbon core, where there is an increased 
availability of 12 C fuel from mixing. During carbon burning, 
the timescale for burning is governed primarily by the neu- 
trino losses (as is true for all advanced burning stages) and 
these losses generally increase monotonically with increas- 
ing temperature. In fact, massive star evolution during the 
advanced stages of evolution can be described as a neutrino- 



mediated Kelvin-Helmholtz contraction of a carbon-oxygen 
core ( jWoosley et~aLl[2002l |E1 Eid et aT]|2l)09| ). Therefore, 
a reduction in the neutrino-losses has the consequence of 
increasing the lifetime of carbon-burning stages. Only the 
carbon shells in models 32CU and 60CU do not show this 
behaviour (see Fig. [7] and [8|. This can be explained by the 
presence of a previous convective carbon core in those mod- 
els, which reduces the abundance of carbon fuel available 
for burning in these shells. Systematic trends during shell 
burning are less clear because of the rather complicated evo- 
lution of the shell structure, but convective shells often form 
at lower temperatures in CI and CU models (see column 7 
in Table [2J, similar to the situation in the core. For carbon 
core burning, on the other hand, there is a clear increase in 
the lifetime with increasing rate, which is shown in Fig. |11| 

The main neutrino processes during carbon burning are 
those caused by pair production and photoneutrino interac- 
tions (IWoosley et al.||2002| lltoh et al. II 19961). It is worth 



noting that the decrease in temperature in the CI and CU 
models is responsible for a larger proportion of neutrinos 
formed by the photoneutrino process rather than pair pro- 
duction. This trend at larger carbon-burning rates is op- 
posite to the trend with initial mass, which favours higher 
temperatures and production of neutrinos by pair produc- 
tion with increasing initial mass. 
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Table 3. Energy generation and neutrino parameters during core carbon burning. For each model the central values of temperature, 
T, density, p, energy generation rates for nuclear burning, £ nuc , and neutrino losses, e„, and percentage fractions of the total neutrinos 
formed by pair production (/ pa ir)i photoneutrino interactions (/phot) an d other processes (/rest) are specified. These parameters are 
determined at the time when the mass fraction of 12 C is half of the value just prior to carbon burning. 
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0.031 


20CU 


0.588 


2.477 x 


10 4 


3.727 x 


10 6 


-2.643 x 


10 s 


41.935 


57.943 


0.122 


25ST 


0.859 


1.439 x 


10 s 


5.176 x 


10 7 


-4.435 x 


10 7 


95.061 


4.917 


0.022 


25CI 


0.690 


3.942 x 


10 4 


2.603 x 


10 8 


-2.975 x 


10 e 


83.475 


16.490 


0.035 


25CU 


0.603 


1.889 x 


10 4 


4.975 x 


10 8 


-4.533 x 


10 s 


58.913 


41.026 


0.061 


32ST 


0.904 


1.313 x 


10 5 


1.360 x 


10 8 


-1.234 x 


10 8 


97.310 


2.680 


0.010 


32CI 


0.711 


3.532 x 


If) 4 


3.682 x 


10 8 


-5.995 x 


10 6 


89.439 


10.543 


0.018 


32CU 


0.621 


1.510 x 


If) 4 


5.725 x 


10 (i 


-9.148 x 


10 s 


74.347 


25.625 


0.028 


60ST 


0.919 


1.106 x 


10 5 


1.900 x 


10 s 


-1.954 x 


10 8 


98.053 


1.941 


0.006 


60CI 


0.725 


3.260 x 


10 4 


5.863 x 


10 8 


-9.442 x 


10 6 


92.247 


7.741 


0.012 


60CU 


0.625 


1.375 x 


10 4 


6.244 x 


10° 


-1.096 x 


10 6 


77.670 


22.309 


0.021 



E 3.5 




Carbon burning rate 

Figure 11. Carbon core burning lifetimes for all models. Note 
that for some models, the core carbon burning is radiative rather 
than convective. The carbon burning lifetime is defined as the 
time for the mass fraction abundance of 12 C to reduce from 0.3 
per cent of its maximum value to a value of 10 — 3 . 



These effects on the central evolution are responsible for 
the different tracks exhibited by the Cf and CU models with 
respect to the ST models in Figures [1] and [5] For the 15 and 
20Mq models the larger cores cause the CI and CU tracks 
to tend towards the higher temperature, lower density side 
of the ST track, but only for the duration the convective core 
is present. When the star moves onto carbon shell burning, 
the core cools and the track returns to the standard curve. 

As explained above, the overlap exhibited by convective 
shells over the ashes of convective carbon cores is due to the 
ignition of carbon that represents the unburnt remainder 
from carbon-core burning. The presence of this remainder 
is caused by the gradual shrinking of the carbon-core near 
the end of the burning stage. This occurs in model 20CI 



and all CU models, except model 15CU where the shell is 
located at the top of the previous convective carbon core. 
The convective carbon shell in the 20CU model (see Fig. 
[6|, however, shows an interesting structure. In this case a 
carbon shell is ignited at a position that overlaps with the 
core and then shortly after an additional shell is ignited 
at the point corresponding to the top of the previous core. 
Because of the unusual structure, the lifetime given in Table 
[2] for the 20CU model, shell 1, is denned from the onset of 
convection to the time it shrinks back up into the second 
shell. 

The presence of overlap with a carbon core has a sig- 
nificant impact on the composition of the shell at the onset 
of convection. Indeed, carbon-core burning ashes, including 
s-process nuclides, will mix out to a position above the rem- 
nant mass and be present in the supernova ejecta. As men- 
tioned above, overlapping shells have previously been noted 
in the literature, but the consequences of overlapping shells 
of this nature are not well studied. The nucleosynthetic con- 
sequences of overlap will be discussed in !j4] 



3.3 Advanced stages beyond carbon burning 

Despite the changes to the stellar structure during carbon 
burning, the evolution of the advanced burning stages in the 
core following carbon-burning seems only slightly affected in 
terms of the convection zone structure, as seen in Fig. [6] [7] 
and [8] but exhibit burning stages with different lifetimes. 
The burning lifetimes for the hydrostatic burning stages are 
presented in Table |4j which are defined for each stage as 
the difference in age from the point where the principal fuel 
for that stage for hydrogen burning, 4 He for helium 
burning, etc.) is depleted by 0.3 per cent from its maximum 
value to the age where the mass fraction of that fuel de- 
pletes below a value of 10 _J , except for carbon burning and 
neon burning, where this value is 10 -3 , and oxygen burn- 
ing, where this value is 10 -2 . These criteria are necessary 
to ensure a lifetime is calculated in those cases where resid- 
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ual fuel is unburnt (such as during oxygen burning in the 
15CU model, where the 16 mass fraction abundance that 
remains unburnt following the end of core oxygen burning is 
~ 3.177 x 10" 3 ) and to ensure that the burning stages are 
correctly separated. The lifetime of the advanced stages is 
relatively sensitive to the mass fractions of isotopes defining 
the lifetime. 

Carbon burning lifetimes are longer for the CI and CU 
rates, as explained in §3.2| but lifetimes for the other ad- 
vanced stages do not show a general trend with the carbon 
burning rate. This lack of trend also applies to the central 
properties, as seen in Fig. [4j where the tracks are modi- 
fied by the enhanced rate models but the modifications do 
not follow a general pattern. In fact, there are examples of 
T c — p c tracks, e.g. the 25CI and 25CU models in Fig. [IJ 
where following the deviation caused by carbon ignition the 
track returns to that of the ST rate (especially for the 15, 
20 and 25 Mq models). The main property determining the 
variations in the lifetime is the central temperature, which 
is linked with the neutrino loss rates. 

The last column of Table [1] shows that the total life- 
time of the star increases slightly with an enhanced carbon 
burning rate, because of the longer carbon burning lifetime. 
Since the total lifetime increases by « 1 — 5 x 10 4 years, 
the strong mass loss (characteristic of massive stars), which 
can increase by up to ~ 10 _5 Mq yr _1 , increases the mass 
lost by up to 0.5Mq. This is demonstrated in column 2 of 
Table [5j which shows the core masses at the end of oxy- 
gen burning for all models. In column 3 of Table [5] we see 
that the carbon burning rate does not affect the helium core 
mass (the helium core mass is defined as the mass coordi- 
nate where the mass fraction abundance of 4 He is 0.75 at 
the interface between the hydrogen and helium-rich layers). 
There is only a tiny difference for the 25 Mq case because of 
the small structure re-arrangement of the hydrogen burning 
shell. In column 4, we see that with an increasing carbon 
burning rate, the CO core mass is larger (the CO core mass 
is defined as the mass coordinate where the 4 He mass frac- 
tion abundance is 10 -3 ). The reason is the following. With 
an increased rate, carbon burning occurs at lower tempera- 
tures where the energy production dominates over neutrino 
cooling and this leads to a stronger carbon core burning in 
a larger convective zone. Thus the carbon burning core pro- 
duces more energy and this leads to a less energetic helium- 
burning shell that is radiative rather than convective, which 
is the case for the ST models. When the He-shell is radiative 
the burning front depletes completely the helium available 
at one mass coordinate and then moves upwards leading to 
a more massive CO core whereas with a convective He-shell, 
the bottom of the shell stays at the same mass coordinate 
since the helium in the convective shell is never completely 
exhausted due to mixing. Note also that the 32 and 60 Mq 
models do not exhibit a value for M™ % . This is because the 
mass loss is strong enough in these WR stars to expel the 
majority of their helium-rich envelopes and the 4 He abun- 
dance is not high enough to satisfy the criterion for M^ 5% . 
In these cases, the helium core mass is taken as the final 
mass, MFi na i (see column 2 of Table |5J. 

As mentioned above, the size of the convective cores 
during neon, oxygen and silicon burnings is only slightly 
affected by the changes in carbon burning rate, as can be 
seen in the last column of Table [5] for the oxygen-free core, 



Mo-free, calculated at the end of core oxygen burning. The 
changes in Mo-frcc with carbon burning rate are because of 
changes in the position of the lower boundary of the last con- 
vective carbon shell. Generally, the magnitude of the changes 
in Mo -free are small and do not present a clear pattern. 



4 NUCLEOSYNTHESIS 
4.1 Neutron sources 

The main effects on the nucleosynthesis in the stellar mod- 
els are due to the lower central temperature of the star 
and the increased lifetime. In particular, the lower cen- 
tral temperature will affect the efficiency of neutron source 
reactions. We recall that the main neutron sources for 
the s process are 13 C, which is important during carbon 
core burning, and 22 Ne, which is important during helium 
core burning and carbon shell burning. The 13 C neutron 
source is mainly produced during carbon core burning by 
the 12 C(p,7) 13 N(/3 + ) 13 C reaction chain. Neutrons are then 
produced by 13 C(a,n) 16 reactions. The protons and a- 
particles originate directly from the 12 C + 12 C fusion re- 
actions. There is competition between the 13 N(/3 + ) 13 C and 
13 N(7,p) 12 C, where at temperatures above 0.8 GK, the (7,p) 
reaction dominates over the /3-decay. The 13 C neutron source 
is thus an efficient neutron producer only at lower tempera- 
tures. During carbon shell burning, where the temperatures 
are higher, the 22 Ne source is the dominant neutron source. 
One can therefore expect that as the carbon burning rate is 
increased and the interior temperature is lowered, the effi- 
ciency of the 13 C neutron source will increase. This efficiency 
will also be higher given the increased lifetimes. 

A non-negligible fraction of neutrons are also present 
from the 17 and 21 Ne neutron sources, but these nuclei are 
mainly produced by neutron captures on 16 O and 20 Ne (and 
17 0(a, 7) 21 Ne) and therefore only act as mediators of the 
neutron irradiance. The 25 Mg(«,n) 28 Si and 12 C( 12 C,n) 23 Mg 
neutron sources are marginal for all models considered here, 
despite the increases to the carbon burning rate. We refer to 
Pignatari et al. (2011) for a more detailed discussion about 
the 12 C( 12 C,n) 23 Mg reaction. 



4.2 S-process parameters 

Several indicators for the neutron capture nucleosynthesis 
are considered. The s process is typically characterised by 
the neutron density, n n , the neutron captures per iron seed, 
n c , and the neutron exposure, r n . n c is defined as follows: 



E «( J 4 i -56)(X i -XP) 



(4.1) 



where Xf is the initial mass fraction abundance of 
isotope Xi with atomic mass Ai and Xs6 Fc is the intial 
mass fraction abundance of 56 Fe, which is the dominant 
seed isotope for s-process nucl eosynthesis. r n is defined as 
T n = / vrrindt (Clayton 19681. However, these definitions 
are of limited use in the multi-zone calculations used here. 
The reason for this is that in the multi-zone stellar models, 
convective mixing affects the neutron irradiance experienced 
by a given mass element ( The et al.||2007 l. Stellar matter, 
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Table 4. Lifetimes for all core burning stages in all models (in yrs). Lifetimes are provided for hydrogen burning (th), helium burning 
( T He)> carbon burning (re), neon burning (tn c ), oxygen burning (tq) and silicon burning (rg;). The total lifetime is given by (iTotal)- 



Model TH THc TC TNc TQ TSi T Total 



15ST 


1 


137 


X 


10 7 


1.255 


X 


10 6 


2 


595 


X 


10 3 


1 


253 


1 


233 


1 


685 


X 


10- 


2 


1 


268 


X 


10 7 


loOl 


1 




X 


1U 




X 


10 


1 


1 00 


X 


-1 r\4 
1U 


11 


one 


1 


1 A ^ 












1 




X 


1 n? 
1U 


15CU 


1 


137 


X 


10 7 


1.255 


X 


10 6 


5 


288 


X 


10 4 


12 


918 


8 


815 












1 


272 


X 


10 7 


20ST 


7 


926 


X 


10" 


8.396 


X 


10 5 


7 


409 


X 


10 2 





193 





293 


1 


302 


X 


10- 


2 


8 


799 


X 


10 6 


20CI 


7 


926 


X 


10 6 


8.396 


X 


10 s 


6 


786 


X 


10 3 





655 





542 












8 


803 


X 


10" 


20CU 


7 


926 


X 


10 s 


8.396 


X 


10 s 


3 


275 


X 


10 4 





265 





253 












8 


825 


X 


10 e 


25ST 


6 


492 


X 


10 6 


6.519 


X 


10 s 


3 


131 


X 


10- 





634 





603 


4 


322 


X 


10" 


3 


7 


168 


X 


10 6 


25CI 


6 


492 


X 


10" 


6.519 


X 


10 5 


2 


984 


X 


10 3 





539 





597 


1 


097 


X 


10- 


2 


7 


169 


X 


10 6 


25CU 


6 


492 


X 


10° 


6.519 


X 


10 s 


2 


296 


X 


10 4 





505 





515 


1 


746 


X 


10- 


2 


7 


186 


X 


10 6 


32ST 


5 


287 


X 


10 8 


5.346 


X 


10 s 


1 


245 


X 


10 2 





111 





167 


8 


997 


X 


10- 


3 


5 


840 


X 


10 6 


32CI 


5 


287 


X 


10 6 


5.346 


X 


10 s 


1 


406 


X 


10 3 





726 


1 


123 


1 


173 


X 


10- 


2 


5 


840 


X 


10 e 


32CU 


5 


287 


X 


10° 


5.346 


X 


10 6 


1 


419 


X 


10 4 





148 





111 


5 


458 


X 


10- 


3 


5 


852 


X 


10 e 


60ST 


3 


549 


X 


10° 


3.935 


X 


10 s 


7 


808 


X 


10 1 





090 





119 


8 


624 


X 


10- 


3 


3 


955 


X 


10" 


60CI 


3 


549 


X 


10 6 


3.935 


X 


10 s 


1 


132 


X 


10 3 





425 





505 












3 


955 


X 


10 6 


60CU 


3 


549 


X 


10 8 


3.935 


X 


10 s 


1 


331 


X 


10 4 





112 





071 












3 


966 


X 


10 e 



Table 5. Core masses at the end of oxygen burning, in solar 
masses. For each model, the final total mass (Mp; na i), helium 
core mass (Af 75< ^), CO core mass (Mqq) and the oxygen- free 
core mass (Mo_f rcc ) are specified. Note that the 32 and 60 Mq 
models expel most of their helium-rich envelopes, consequently 
becoming WR stars. 



Model 




M 75% 




A^O-trco 


15ST 


12.132 


4.791 


2.805 


0.921 


15CI 


12.069 


4.791 


2.923 


0.867 


15CU 


11.907 


4.791 


3.239 


0.849 


20ST 


13.974 


6.826 


4.494 


1.083 


20CI 


13.916 


6.826 


4.491 


1.099 


20CU 


13.602 


6.826 


4.696 


1.040 


25ST 


13.738 


9.199 


6.301 


1.081 


25CI 


13.710 


9.092 


6.384 


0.980 


25CU 


13.202 


9.092 


6.544 


1.124 


32ST 


12.495 


12.495 


9.146 


1.187 


32CI 


12.495 


12.495 


9.146 


0.984 


32CU 


12.493 


12.493 


9.425 


1.334 


60ST 


13.428 


13.428 


10.701 


1.242 


60CI 


13.423 


13.423 


10.446 


0.990 


60CU 


13.278 


13.278 


10.929 


1.519 



including the neutron sources, seeds and poisons, is mixed 
into and out of the bottom of the convection zone, where 
the temperature is highest and where the majority of the s 
process occurs. Consequently, an evaluation of n c or r n at a 
particular mass coordinate will be different to that experi- 
enced by a given mass element. 

Therefore, in order to evaluate relevant parameters to 
describe the neutron irradiance, convective mixing needs to 
be taken into account in the evaluation of the parameter. 
This can be achieved for the neutron exposure by consider- 
ing the initial and final abundances of 54 Fe, an isotope that is 
slowly destroyed by neutron captures in the s-process sites 
considered here. It cannot be used during or after oxygen 
burning where temperatures are high enough to photodisin- 
tegrate heavy elements ( Woosley &: Weaver||1995| ). An esti- 
mate of the neutron exposure using M Fe can be made using 



the following formula (Woosley & Weaver 1995 
2000|, 



The et al. 



T54 



; [lnXi( 54 Fe) -lnX f ( 54 Fe)], 



(4.2) 



where 
mb, 



a is the 54 Fe(n,7) 55 Fe reaction rate (a = 29.6 ± 1.3 
, Dillmann et~al7||2006| ) and X ; ( 54 Fe) and X f ( 54 Fe) are 
the mass fraction abundances of 54 Fe before and after the 
neutron exposure respectively. A better estimate of n c can 
be obtained by using mass-averaged abundances for X iy Xf 
and X56 Fe over the maximum size of the convective region, 



Er(^-56)((X i )-(X I °)) 

(A56 Fc ) 



(4.3) 



This takes into account any changes to the size of the con- 
vective region during the burning stage where the s-process 
nucleosynthesis occurs. 

Table [5] lists, for all models, the neutron exposure, T54, 
the neutron captures per iron seed, n Ciav , the mass fraction 
abundances of the isotopes 54 Fe and 88 Sr and the isobaric 
ratios 70 Ge/ 70 Zn, 80 Kr/ 80 Se and 86 Sr/ 86 Kr. 88 Sr, like 54 Fe, 
is also a useful s-process indicator as it has a neutron-magic 
nucleus (iV = 50) and is slowly built-up over the course of 
the s process. The isobaric ratios are also specified, because 
changes to the ratios demonstrate deviations to the s-process 
path at branching point nuclides ( 69 Zn, 79 Se and 85 Kr for 
70 Ge/ 70 Zn, 80 Kr/ 80 Se and 86 Sr/ 86 Kr respectively). Indeed, 
if the neutron density increases, the s-process path opens to 
allow the production of more neutron-rich isotopes, lowering 
these ratios. 



4.3 Core carbon burning 

According to Table [6] all CI and CU models show a de- 
pletion of o4 Fc and production of 88 Sr relative to the ST 
case, indicating that a higher neutron exposure is present 
in the convective carbon core. For all CI and CU models, 
irrespective of mass, the neutron exposure is high enough 
to allow an increasing production of isotopes beyond the 
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Table 6. S-process tracers, neutron capture parameters and isotopic ratios at the end of helium-core burning, carbon-core burning 
and convective carbon-shell burning. n c , a v is the neutron captures per iron seed averaged over the convective region and T54 is the 
neutron exposure calculated using Eq. 4.2 The 88 Sr and 54 Fe abundances are specified as average mass fraction abundances, -X"s8g r and 



^54 Pc respectively, at the end of the burning stage over the convective region, except for radiative burning where the central values are 
taken. The s-process parameters for a shell that persists to the presupernova stage use final abundances that are evaluated at start of 
oxygen burning, which removes the effects of photodisintcgration occurring during the late evolutionary stages from the evaluation of 
the s-process parameters. 



Model 


Shell 


8 


3 Sr 




54 Fe 






TKA ( TTlV> ^ 


70 Ge/ 70 Zn 


80 Kr/ 80 Se 


86 Sr/ 86 Kr 


15ST 


He-core 


2.005 


X 


10" 


7 


5.750 x 


10 


-6 


1.641 


0.088 


115.913 


2.690 


4.247 


15ST 


C-core 


1.556 


X 


10- 


-6 


7.721 x 


10 


-7 


6.601 


0.062 


1165.633 


5.107 


46.001 


15ST 


1 


1.000 


X 


10" 


-6 


1.089 x 


10 


-6 


4.740 


0.048 


1036.915 


3.668 


20.178 


15ST 


2 


6.629 


X 


nr 


7 


1.266 x 


10 


-6 


3.903 


0.042 


335.818 


0.701 


2.708 


15CI 


C-core 


1.009 


X 


10- 


4 


9.137 x 


10 


-8 


29.270 


0.134 


901.882 


4.284 


45.048 


15CI 


1 


2.803 


X 


10- 


-5 


6.958 x 


10 


-7 


6.005 


0.059 


862.687 


3.172 


23.268 


15CU 


C-core 


2.182 


X 


10- 


-4 


3.716 x 


10 


-8 


46.293 


0.165 


743.822 


4.080 


44.065 


15CU 


1 


5.046 


X 


10- 


5 


2.163 x 


10 


-6 


19.423 


0.055 


638.189 


0.765 


1.726 


20ST 


He-core 


3.817 


X 


10- 


7 


1.070 x 


10 


-6 


3.069 


0.143 


928.859 


3.588 


7.503 


20ST 


C-core 


1.286 


X 


10- 


-6 


1.615 x 


10 


-7 


8.080 


0.062 


1315.250 


4.012 


30.741 


20ST 


1 


1.064 


X 


10" 


-6 


2.303 x 


10 


-7 


6.605 


0.043 


1245.114 


2.605 


17.583 


20ST 


2 


9.403 


X 


10- 


7 


3.382 x 


10 


-7 


4.934 


0.033 


518.314 


0.774 


4.205 


20ST 


3 


8.762 


X 


10- 


7 


4.292 x 


10 


-7 


0.119 


0.001 


487.403 


0.696 


3.802 


20CI 


C-core 


5.197 


X 


10- 


-5 


8.818 x 


10 


-8 


27.796 


0.084 


970.039 


4.200 


41.853 


20CI 


1 


2.424 


X 


10- 


-5 


3.828 x 


10 


-7 


5.920 


0.023 


975.182 


2.873 


20.450 


20CI 


2 


2.160 


X 


10- 


5 


3.869 x 


10 


-7 


2.737 


0.012 


347.183 


0.366 


3.352 


20CU 


C-core 


1.727 


X 


10- 


4 


4.802 x 


10 


-9 


60.722 


0.182 


779.749 


4.104 


36.648 


20CU 


1 


7.074 


X 


10- 


-5 


5.484 x 


10 


-7 


4.073 


0.019 


494.139 


2.019 


22.567 


20CU 


2 


1.194 


X 


10- 


-5 


6.573 x 


10 


-7 


4.651 


0.027 


151.579 


0.348 


4.048 


25ST 


He-core 


6.153 


X 


10- 


7 


3.539 x 


10 


-7 


4.280 


0.180 


2220.036 


3.755 


11.329 


25ST 


C-core 


1.472 


X 


10- 


6 


7.918 x 


10 


-8 


8.271 


0.045 


1432.597 


4.385 


35.554 


25ST 


1 


9.499 


X 


10- 


-7 


1.482 x 


10 


-7 


5.632 


0.028 


87.609 


0.109 


0.515 


25CI 


C-core 


4.092 


X 


10- 


-5 


1.411 x 


10 


-9 


48.421 


0.179 


970.416 


4.576 


59.426 


25CI 


1 


1.772 


X 


10- 


5 


6.313 x 


10 


-8 


23.538 


0.045 


1063.729 


4.066 


38.990 


25CI 


2 


1.111 


X 


10- 


6 


1.564 x 


10 


-7 


5.543 


0.028 


315.357 


0.280 


1.401 


25CU 


C-core 


1.475 


X 


10- 


-4 


1.509 x 


10 


-9 


73.339 


0.184 


804.018 


4.072 


36.419 


25CU 


1 


9.824 


X 


10- 


-5 


1.347 x 


10 


-7 


15.755 


0.015 


698.157 


1.283 


10.094 


32ST 


He-core 


1.097 


X 


10- 


6 


1.192 x 


10 


-7 


5.623 


0.217 


3380.614 


3.900 


16.340 


32ST 


C-core 


1.788 


X 


10- 


6 


5.333 x 


10 


-8 


6.239 


0.024 


1640.445 


3.640 


28.449 


32ST 


1 


1.315 


X 


10- 


-6 


8.625 x 


10 


-8 


3.016 


0.010 


75.996 


0.130 


1.014 


32CI 


C-core 


1.825 


X 


10- 


-5 


3.955 x 


10 


-9 


38.296 


0.110 


1042.993 


4.740 


60.126 


32CI 


1 


2.045 


X 


10- 


6 


6.562 x 


10 


-8 


5.220 


0.017 


1021.836 


1.646 


9.944 


32CU 


C-core 


1.007 


X 


10- 


4 


8.498 x 10- 


10 


77.718 


0.167 


837.791 


3.949 


39.032 


32CU 


1 


7.633 


X 


10- 


-5 


3.346 x 10 


-8 


16.738 


0.011 


509.651 


0.428 


4.911 


60ST 


He-core 


1.524 


X 


10- 


-6 


6.404 x 10 


-8 


6.489 


0.238 


1741.270 


1.125 


12.267 


60ST 


C-core 


1.701 


X 


10- 


6 


5.297 x 10 


-8 


5.862 


0.023 


1743.568 


3.246 


25.865 


60ST 


1 


1.335 


X 


10- 


6 


7.814 x 10 


-8 


2.779 


0.009 


69.670 


0.146 


1.136 


60CI 


C-core 


1.491 


X 


10- 


-5 


4.808 x 10 


-9 


33.897 


0.104 


1072.384 


4.619 


52.637 


60CI 


1 


1.622 


X 


10- 


-6 


5.837 x 10 


-7 


3.800 


0.029 


871.777 


0.921 


5.676 


60CU 


C-core 


1.076 


X 


10- 


4 


6.551 x 10- 


10 


81.743 


0.172 


837.512 


3.877 


36.865 


60CU 


1 


8.908 


X 


10- 


5 


2.512 x 10 


-8 


17.940 


0.010 


455.999 


0.370 


4.862 



Sr-Y-Zr peak, which is quantified in a higher neutron cap- 
tures per iron seed. An example of this nucleosynthesis for 
the 15Mq model is seen in Fig. |12[ which shows the central 
overproduction factors for heavy, stable isotopes in the star 
at the end of carbon burning. The distribution of synthe- 
sised isotopes is extended, with increasing rate, beyond the 
Sr-Y-Zr peak to include isotopes up to the Ba-La peak at 
A ~ 140. This is an anomalous distribution compared to the 
weak s-process component. 

The neutron density in the carbon core decreases from 
a typical value of ~ 10 8 cm -3 , which is maintained through- 
out the burning, to ~ 10 7 cm~ 3 in the models with an in- 
creasing carbon burning rate. In the 25CU, 32CU and 60CU 



models the neutron density is enhanced over the CI cases be- 
cause of the presence of the convective core; the mixing into 
and out of the centre acts to maintain a supply of neutron 
sources at the centre. Concerning the ST case, the neutron 
exposures for the cores are similar in magnitude to that of 
the helium burning core (~ 0.06 mb _1 ), but are lower for 
the most massive stars considered here (~ 0.02 mb _1 for 
the 32ST and 60ST models). For the CI and CU rates the 
neutron exposures are significantly enhanced, typically ex- 
ceeding 0.1 mb _1 . This is mainly due to the rising efficiency 
of the 13 C neutron source at lower temperatures, coupled 
with the increased lifetime of the core carbon burning stage. 
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Figure 12. Central overproduction factors for most stable isotopes at the end of central carbon burning for the 15Mq models. The plot 
shows a significant increase in nucleosynthesis of isotopes between 60 < A < 140 in the CI and CU models, which is beyond the Sr-Y-Zr 
peak at an atomic mass RJ 90. 



4.4 Carbon shell burning 

Nucleosynthesis in the carbon shells is characterised by the s 
process with a high neutron density but lower neutron expo- 
sure compared to carbon core, with 22 Ne being the dominant 
neutron source. In the ST models, the neutron densities vary 
from ~ 10 s cm -3 for early convective shells (models 15ST 
and 20ST), and increase to a typical value of ~ 10 10 cm -3 in 
the final carbon burning shell. In the CI and CU models, the 
neutron density is ~ 10 7 cm -3 in early shells, similar to the 
values obtained during core carbon burning, and then rises 
to ~ 10 8-9 cm -3 . The lifetimes for the carbon shell burning 
stages vary quite differently from model to model, but are 
generally increasing with increasing rate. For example, in 
the 15CU case, the lifetimes of the last carbon shell in Ta- 
ble[2]for the 15ST, CI and CU models are 17.92, 150.1 and 
594.2 years respectively. The carbon shell in model 15CU 
consequently exhibits a strong neutron exposure of similar 
magnitude to the carbon core (see Table |6|. It should be 
noted however that in almost every instance of a carbon 
burning shell, the neutron exposure is smaller than that of 
the carbon core in the same model. This asserts the fact 
that carbon shells are characterised by a lower neutron ex- 
posure and higher neutron density (with 22 Ne as the main 
neutron source), although the degree with which this is true 
is reduced with an increasing carbon burning rate. That is, 
the general trend with increasing rate is a decrease in the 
neutron density and an increase in the neutron exposure in 
the carbon shells. 

The above can be verified by considering the ratios of 
isotopes involved at branching points, since the lower neu- 
tron density will close the s-process path to the synthesis 



of more neutron-rich isotopes at branching points. The last 
three columns of Table [5] show the isobaric ratios at the end 
of the core and shell carbon burning stages for 70 Ge/ 70 Zn, 
80 Kr/ 80 Se and 86 Sr/ 86 Kr, with values for the end of helium 
core burning specified for reference. For most models, the 
ratios increase in the last carbon shell with increasing car- 
bon burning rate, favouring production of the s-only isotopes 
70 Ge, 80 Kr and 86 Sr, due to the lower neutron density in the 
carbon-shells in the CI and CU models. However, the ratios 
are sensitive to convection, since shell overlap causes the 
shells to be polluted with carbon core s-process ashes. Con- 
sequently, the 25CU, 32CU and 60CU models instead show 
a decrease in the ratios. Considering that the ratios in the 
initial composition are 3.271, 6.124 and 0.036 for 70 Ge/ 70 Zn, 
80 Kr/ 80 Se and 86 Sr/ 86 Kr, the presence of lower isobaric ra- 
tios than these in the shells indicates that the branching is 
indeed affected during the carbon shell s-process and that 
the decrease is not associated purely with the mixing of car- 
bon core matter with helium burning ashes. 



5 YIELDS 

5.1 Calculations 

The yields calculations were made in the same manner as 
that of Hirschi et al. (20051, which considers two contri- 



butions to the yields: the stellar wind and the supernova 
explosion. The wind yield for nuclide i for a star with initial 
mass m is calculated using: 
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wind 

mp im 



r(m) 



M(m,t)[X?(m,t) - 



X°]dt 



(5.1) 



where r(m) is the final age of the star, M(m, t) is the mass 
loss rate, Xf is the surface mass-fraction abundance, Xf 
is the initial mass-fraction abundance. The majority of the 
matter lost through the stellar wind occurs during hydrogen 
and helium burning. The composition of the wind is similar 
to that of the initial composition, except for the 32Mq and 
60Mq models where the mass loss is significant enough to 
include some of the hydrogen burning ashes. Table [5] shows 
that the total mass lost over the stellar evolution due to the 
stellar wind increases significantly with initial mass (~ 20 
per cent lost for the 15 Mq models to « 80 per cent lost for 
the 60 Mq models). 

The presupernova yields are calculated using: 



prcSN 
m Pim 



[Xi(m r ) — Xi]dm r 



(5.2) 



where m T is the total mass of the star at r(m), M rem , m is 
the remnant mass, Xf is the initial mass fraction abundance 
of element i and Xi(m r ) is the mass fraction abundance at 
mass coordinate m r . The total yields are then just the sum of 
the wind and the presupernova yields. The calculated yields 
of selected isotopes for model 15ST are shown in Table[7] (full 
yield tables for all models are provided with the electronic 
edition of this paper). 

The point in the evolution in which the yields are taken 
in this work is at the end of central oxygen burning, as ex- 
plained in §2.2| This choice was made since not all the mod- 
els were post-processed until the end of silicon burning. No- 
tice that, as mentioned in §2.2| after central oxygen burning, 
the material outside the remnant mass is not affected much 
by the pre-explosive evolution. The only potential contribu- 
tions that may affect the s-process abundances are during 
the early collapse, when the neutron density may increase 



significantly (e.g. in the carbon shell, see Pignatari et al 
2010[ ) or partial or complete photodisintegration at the bot- 
tom of the carbon, neon and oxygen shells. The effects of 
photodisintegration will be discussed in a forthcoming pa- 
per (Pignatari et al. 2011, in prep.). 

With regards to explosive burning, the supernova explo- 
sion is responsible for destroying and recreating a portion of 
the ejecta, which includes p-process rich and, to a smaller 
extent, s-process rich layers, possibly having a relevant im- 
pact on the total yields of s-process nuclides (see for instance 
Rauscher et al.||2002| |Tur et al.|[2009[ ). However, the explo- 
sive burning process is sensitive to uncertainties in the su- 
pernova explosion mechanism for the range of initial masses 



considered here (Fryer 20091. The uncertainties associated 



with the supernova explosion, namely the explosion energy, 
the ignition mechanism and the amount of fall-back, are im- 
portant especially for the 15, 20 and 25 Mq models. These 
uncertainties would also affect the amount of matter locked 
up in the remnants. In this work, the remnant mass takes 
into account the additional matter that falls back onto the 
remnant following the initial explosion. The choice of rem- 
nant masses for the models is taken from the analytical fits 
of Fryer et al. (2011, in prep.) for solar metallicity stars, 



which derive from energy-driven explosions (see for instance 
Fryer|20 09). The remnant masses, Af rC m,m, are given by 



M 



rcm.m 



1.1 + 0.2e 



(m-ll)/4 



3e 



0.4(m-26) 



18.35 -0.3m, 



11 < m ^ 30 
30 < m < 50 
(5.3) 

which gives remnant masses of 1.61, 2.73, 5.71 and 8.75 Mq 
for initial masses, m, of 15, 20, 25 and 32 Mq respectively. 
For the 60 Mq models a remnant mass was calculated by 
scaling with the CO core mass ratio for the ST models, 



M, 



Mr, 



McO,60Mq 
© I McO,32M, 



© 



(5.4) 



giving a remnant mass of 10.24 Mq. The resultant remnant 
masses are such that for the 15 Mq models, the oxygen 
shell is partially included in the supernova ejecta. For the 
other models however, the remnants are large and the ejecta 
includes the upper portion of the carbon shell and the over- 
lying layers only. The remnant masses here are larger in 
comparison with those used in previous studies of explosive 
nucleosynthesis ( Limongi et al.|2000 Rauscher et al.|2002 |. 
This is due to the use, in those studies, of piston-driven mod- 
els that are known to underestimate the amount of fall-back 
onto the supernova remnant (Young & Fryer 20071. The 



large remnant masses may cause the explosive nucleosyn- 
thesis to occur predominantly in the layers that fall back 
onto the remnant. 

In addition to the yields, the ejected masses, Ei m can be 
calculated, which are the exact analogues of Eq. |5.1| and |5.2| 
but without the inclusion of the Xf term. If the total mass of 
matter ejected is M c - hm = m T — M romjm , the overproduction 
factors averaged over the ejecta are calculated using 



(OP) lm = 



E 

^ in 



MeUmXf 



(5.5) 



The overproduction factors averaged over the ejecta for 
the s-only isotopes are shown in Fig. |13| which represents 
well the general abundance distribution for stable isotopes 
created by the models. A considerable amount of s-process 
nucleosynthesis occurs for all CU models by up to 3 dex, 
which is either because of overlap between the carbon shells 
and the carbon core (for models 20CI, 25CU, 32CU and 
60CU) or because of strong neutron exposures in the carbon 
shells (models 15CU and 20CU). The 20CI model features a 
strong overlap between the convective carbon core and the 
successive carbon shells, which is not seen in model 20CU 
and therefore has more significant production than model 
20CU. In fact, for the CI rate, only the 20 Mq model shows 
a significantly enhanced production over the ST rate. The 
15CI model also shows some production, but the distribution 
of isotopes is very similar to that of model 15ST. This is 
in contrast to the 20CI model, which shows an extended 
distribution of production featuring heavier nuclides. 

A first order approximation of the weak s-process com- 
ponent can be made by taking the sum of the yields for each 
stellar model, taking into account the number of stars with 
that initial mass formed, 



V 



(5.6) 
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Table 7. Yields for model 15ST. For each isotope, i, the atomic mass (A), atomic number (Z), initial mass fraction abundance (Xf), 
wind yield (mp wmd , in Mq), presupernova yield (mp preSN , in Mq), total yield (mp total , in Mq), total ejected mass (Ei m , in Mq) and 
average overproduction factor ({OP)) arc specified. The decays of unstable species to their stable isobars are taken into account. 



Isotope 


A 


Z 


X? 


mp wind 


mp preSN 


m ptotal 


R- 


(OP) 





1 


1 


7 


.064E-01 


-4 


366E-02 


-2.933E+00 


-2.977E+00 


6.485E+00 





.685 


4 He 


4 


2 


2 


.735E-01 


4 


345E-02 


1.435E+00 


1.479E+00 


5.142E+00 


1 


.404 


12 C 


12 


6 


3 


.425E-03 


-2 


.639E-03 


3 


101E-01 


3 


.074E-01 


3 


.533E-01 


7 


.703 


13 C 


13 


6 


4 


.156E-05 


2 


.302E-04 


2 


276E-04 


4 


.577E-04 


1 


.014E-03 


1 


.822 


14 N 


14 


7 


1 


.059E-03 


4 


.132E-03 


3 


.401E-02 


3 


.814E-02 


5 


.232E-02 


3 


.689 


16Q 


16 


8 


9 


.624E-03 


-1 


474E-03 


7. 


579E-01 


7 


.564E-01 


8 


.853E-01 


6 


.868 


19p 


19 


9 


5 


.611E-07 


-9 


796E-08 


-2. 


190E-06 


-2 


.288E-06 


5 


.227E-06 





.696 


20 Ne 


20 


10 


1 


.818E-03 


-2 


.514E-06 


3 


238E-01 


3 


.238E-01 


3 


.482E-01 


14 


.302 


23 Na 


23 


11 


4 


.000E-05 


3 


.023E-05 


1 


337E-02 


1 


.340E-02 


1 


.394E-02 


26 


.021 


24 Mg 


24 


12 


5 


.862E-04 


-1 


.079E-08 


2. 


.747E-02 


2 


.747E-02 


3 


.532E-02 


4 


.498 


27 A1 


27 


13 


6 


.481E-05 


4 


579E-08 


3 


142E-03 


3 


.142E-03 


4 


.010E-03 


4 


.620 


28 Si 


28 


14 


7 


.453E-04 


-1 


.752E-08 


1 


.844E-03 


1 


.844E-03 


1 


.183E-02 


1 


.185 


31p 


31 


15 


7 


.106E-06 


1 


.394E-09 


7. 


106E-05 


7 


.106E-05 


1 


.662E-04 


1 


.747 


32g 


32 


16 


4 


.011E-04 


-9 


512E-09 


-1 


897E-04 


-1 


.897E-04 


5 


.182E-03 





.965 


36 Ar 


36 


18 


8 


.202E-05 


-1 


.944E-09 


-7. 


472E-05 


-7 


.472E-05 


1 


.024E-03 





.932 


39 K 


39 


19 


3 


.900E-06 


-9 


.244E-11 


7. 


.466E-06 


7 


.466E-06 


5 


.970E-05 


1 


.143 


40 Ca 


40 


20 


7 


.225E-05 


-1 


.706E-09 


-5 


212E-05 


-5 


.212E-05 


9 


.156E-04 





.946 


45 Sc 


45 


21 


5 


.414E-08 


-1 


.283E-12 


8 


303E-07 


8 


.303E-07 


1 


.555E-06 


2 


.145 


50 T ; 


50 


22 


2 


.208E-07 


-5 


.234E-12 


3 


801E-06 


3 


.801E-06 


6 


.758E-06 


2 


.285 


Sly 


51 


23 


4 


.138E-07 


-9 


.808E-12 


-6. 


535E-08 


-6 


.536E-08 


5 


.476E-06 





.988 


52 Cr 


52 


24 


1 


658E-05 


-3 


.929E-10 


-1 


.282E-05 


-1 


.282E-05 


2 


.092E-04 





.942 


55 Mn 


55 


25 


1 


.098E-05 


-2 


603E-10 


3. 


666E-06 


3 


.666E-06 


1 


.507E-04 


1 


.025 


54 Fe 


54 


26 


8 


.118E-05 


-1 


.924E-09 


-1 


.208E-04 


-1 


.208E-04 


9 


.665E-04 





.889 


56 Fe 


56 


26 


1 


.322E-03 


-3 


.133E-08 


-1 


213E-03 


-1 


.213E-03 


1 


.649E-02 





.931 


59 Co 


59 


27 


3 


991E-06 


-9 


461E-11 


2 


.580E-04 


2 


.580E-04 


3 


.114E-04 


5 


.825 


60 Ni 


60 


28 


2 


.276E-05 


-5 


.394E-10 


1. 


437E-04 


1 


.437E-04 


4 


.485E-04 


1 


.472 


63 Cu 


63 


29 


6 


.600E-07 


-1 


.564E-11 


5 


.493E-05 


5 


.493E-05 


6 


.376E-05 


7 


,213 


65 Cu 


65 


29 


3 


.035E-07 


-7 


.193E-12 


3 


249E-05 


3 


.249E-05 


3 


.655E-05 


8 


.993 


64 Zn 


64 


30 


1 


.131E-06 


-2 


.680E-11 


1 


.792E-05 


1 


.792E-05 


3 


.306E-05 


2 


183 


66 Zn 


66 


30 


6 


.690E-07 


-1 


586E-11 


1 


.856E-05 


1 


.856E-05 


2 


.752E-05 


3 


.072 


70 Zn 


70 


30 


1 


.577E-08 


-3 


737E-13 


-1. 


160E-08 


-1 


.160E-08 


1 


.996E-07 





.945 


69 Ga 


69 


31 


4 


.551E-08 


-1 


.079E-12 


2 


367E-06 


2 


.367E-06 


2 


.977E-06 


4 


.884 


71 Ga 


71 


31 


3 


.108E-08 


-7 


.366E-13 


2 


012E-06 


2 


.012E-06 


2 


.428E-06 


5 


.834 


70 Ge 


70 


32 


5 


.157E-08 


-1 


.222E-12 


3 


.185E-06 


3 


.185E-06 


3 


.876E-06 


5 


.611 


72 Ge 


72 


32 


6 


.910E-08 


-1 


638E-12 


2 


.614E-06 


2 


.614E-06 


3 


.539E-06 


3 


.824 


75 As 


75 


33 


1 


.430E-08 


-3 


.390E-13 


4. 


113E-07 


4 


.113E-07 


6 


.028E-07 


3 


.147 


76 Se 


76 


34 


1 


.296E-08 


-3 


072E-13 


6 


.260E-07 


6 


.260E-07 


7 


.995E-07 


4 


.606 


78 Se 


78 


34 


3 


.376E-08 


-8 


003E-13 


1 


.441E-06 


1 


.441E-06 


1 


.894E-06 


4 


.188 


80 Se 


80 


34 


7 


.226E-08 


-1 


713E-12 


2 


.985E-07 


2 


.985E-07 


1 


.266E-06 


1 


.308 


79 Br 


79 


35 


1 


.389E-08 


-3 


293E-13 


1 


867E-07 


1 


.867E-07 


3 


.728E-07 


2 


.003 


8i Br 


81 


35 


1 


.386E-08 


-3 


.285E-13 


2 


041E-07 


2 


.041E-07 


3 


.897E-07 


2 


.100 


80 Kr 


80 


36 


2 


.575E-09 


-6 


.103E-14 


2 


.610E-07 


2 


.610E-07 


2 


.955E-07 


8 


.569 


82 Kr 


82 


36 


1 


.320E-08 


-3 


.128E-13 


7 


.028E-07 


7 


.028E-07 


8 


.795E-07 


4 


.977 


84 Kr 


84 


36 


6 


.602E-08 


-1 


565E-12 


1 


031E-06 


1 


.031E-06 


1 


.915E-06 


2 


.166 


86 Kr 


86 


36 


2 


.044E-08 


-4 


.846E-13 


1 


289E-07 


1 


.289E-07 


4 


.027E-07 


1 


.471 


85 Rb 


85 


37 


1 


.282E-08 


-3 


040E-13 


1. 


721E-07 


1 


.721E-07 


3 


.438E-07 


2 


.002 


87r>u 
°'Hb 


87 


37 


5 


063E-09 


-2 


025E-12 


6 


776E-08 


6 


.776E-08 


1 


.356E-07 


1 


.999 


84 Sr 


84 


38 


3 


.228E-10 


-7 


651E-15 


-6 


.777E-10 


-6 


.777E-10 


3 


.646E-09 





.843 


86 Sr 


86 


38 


5 


.845E-09 


-1 


.385E-13 


3 


.642E-07 


3 


.642E-07 


4 


.424E-07 


5 


.652 


87 Sr 


87 


38 


4 


.443E-09 


1 


800E-12 


1 


.858E-07 


1 


.858E-07 


2 


.453E-07 


4 


.123 


88 Sr 


88 


38 


5 


011E-08 


-1 


.188E-12 


5 


602E-07 


5 


.602E-07 


1 


.231E-06 


1 


,835 


89y 


89 


39 


1 


.229E-08 


-2 


.914E-13 


9 


875E-08 


9 


.875E-08 


2 


.634E-07 


1 


.600 


90 Zr 


90 


40 


1 


.534E-08 


-3 


637E-13 


4. 


445E-08 


4 


.445E-08 


2 


.500E-07 


1 


.216 


92 Zr 


92 


40 


5 


227E-09 


-1 


.239E-13 


1 


.871E-08 


1 


.871E-08 


8 


.872E-08 


1 


.267 


94 Zr 


94 


40 


5 


.413E-09 


-1 


.283E-13 


6 


178E-09 


6 


.178E-09 


7 


.868E-08 


1 


.085 


93 Nb 


93 


41 


1 


.900E-09 


-4 


504E-14 


7. 


083E-09 


7 


.082E-09 


3 


.253E-08 


1 


.278 


92 Mo 


92 


42 


1 


.012E-09 


-2 


.400E-14 


-1 


.687E-09 


-1 


.687E-09 


1 


.187E-08 





.876 


94 Mo 


94 


42 


6 


.448E-10 


-1 


.528E-14 


2. 


073E-11 


2 


.072E-11 


8 


.656E-09 


1 


.002 


96 Mo 


96 


42 


1 


.188E-09 


-2 


.815E-14 


3 


811E-09 


3 


.811E-09 


1 


.972E-08 


1 


.240 


98 Mo 


98 


42 
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Figure 13. The overproduction factors averaged over the total ejected mass for s-only nuclides as a function of atomic mass. The ST, CI 
and CU rates are indicated by the blue crosses, red diamonds and green circles respectively. Isotopes of the same element are connected 
by adjoining lines. 



where r m is a weighting factor determined by the integra- 
tion of the Salpeter initial mass function (IMF), dN/dm = 
{owi -2,35 , over a certain range. Yields from the 15, 20, 25, 
32 and 60 Mq models were applied to stars within the ini- 
tial mass ranges of 12.5-17.5, 17.5-22.5, 22.5-28.5, 28.5-46 
and 46-80 Mq respectively, giving values of r m equal to 
39.75, 19.89, 13.45, 14.59, 12.32 per cent respectively (with 
£o = 0.304). Consequently, the 15 and 20 Mq models domi- 
nate as the main contributors to the evaluation of the weak 
component (~ 60 per cent of all stars in the total massive 
star mass range considered here). Stars with initial masses 
less than 12.5 Mq or greater than 80 Mq are assumed to 
have a zero contibution to the weak s-process component. 

The 13 C neutron source during carbon core burning is 
mainly primary whereas the 22 Ne source is secondarjQ since 
it depends on the initial 14 N abundance from the CNO cycle. 
If a solar metallicity star of a given mass is the dominant site 
for the production of particular primary and secondary nu- 
clides, A and B, respectively, the overproduction factor for 
B is expected to be approximately twice that of A ( |Truran| 
fc Cameron) 1971 1. Although this is a rather crude approxi- 
mation regarding the detailed nature of chemical evolution 
within galaxies and/or star clusters and the nucleosynthe- 
sis processes themselves (Tinslcy 1979), the weak s process 



5 The products of nucleosynthesis processes in stars, to first or- 
der, can be described as being primary or secondary depending 
on whether the processes responsible for the production depend 
on the initial metallicity. The production of primary nuclides does 
not vary with metallicity whereas secondary nuclides will be pro- 
duced in proportion to their initial seed nuclei. 



in massive stars is expected to hold reasonably to this ap- 
proximation because the dominant neutron sources, seeds 
and poisons of the weak s-process are secondary. It can be 
expected therefore that the overproduction factors for the 
weak s-process nuclides reproduce the solar system abun- 
dances when the overproduction f actor is approximately 
twice that of ie O (Tur et al. 20091. In any case, this rule 



of thumb can be used as a rough guide to indicate the typ- 



ical solar production of s-process nuclides (Rauscher et al. 
2002} |Pignatari et al.|2010 i. 



The overproduction factors of the weak component, 
2/wcak,i/^i°, for nuclides with atomic masses 50 < A < 150, 
are displayed in Fig. |14| Concerning the CU rate, the over- 
production factors are very large (up to 2.56 dex for 86 Sr) 
with respect to the ST model, with significant s-process pro- 
duction of nuclides up to the Ba-La peak at A ss 140. The 
resulting s-process distribution, peaked at the Sr-Y-Zr, is not 
characteristic of the weak s-process component, stopping at 
A « 90. The s-process nuclides with 90 < A < 110 have 
overproduction factors that are comparable to 16 O multi- 
plied by two. Such differences for the CU case compared to 
the classical weak s-process component occur because of the 
13 C neutron source. 

For the CI case, the overabundances of many nuclides 
are similar to the ST case, except for nuclides that are close 
to the Sr-Y-Zr peak or with higher atomic mass (Mo, Ru, 
Cd, and Pd for example). S-process isotopes of Kr and Sr 
have overproduction factors that are higher than 16 O multi- 
plied by two. The abundances of the heavier nuclides Y, Zr, 
Mo, Ru, Cd and Pd show an enhanced production, which is 
0.5 to 1.0 dex lower than the Kr-Sr peak. Overall, the re- 
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suiting s-process distribution is approximately flat from Ni 
to Sr. 

Fig. [15] shows the overproduction factors for the weak 
components of the CI and CU cases plotted relative to the 
ST case. The peak of the relative production of s-process 
nuclides lies at 87 Sr in both cases and declines smoothly 
with increasing mass number, although the overproduction 
factor fo r 86 Sr is slightly larger than the 87 Sr for all cases 
(see Fig. 14 1. For the CU case, the overabundance of 87 Sr is 



1.7 dex larger than for the ST case. The enhancement stops 
at Ba, with 0.5 dex more production and declines steeply, 
with a production of heavier nuclides similar to that of the 
ST case. For the CI case however, the peak production at 
87 Sr is 0.6 dex larger than the ST case and tends to 0.0 at 
Ba. 

The overproduction factors of Sr, Y, Zr, Mo, Ru, Pd 
and Cd are enhanced in the carbon-core s-process (for ex- 



ample, see Fig. 12 1. In the CU case, this occurs for all models 



other than model 15CU. In the CI case, the overlap between 
the convective carbon core and the carbon shell only oc- 
curs for model 20CI. Removing the 20 Mq models from the 
evaluation of the weak component allows for a comparison 
between the predicted weak component with and without 
the occurrence of an overlap. Fig. [16] shows the predicted 
weak component (CI - no20) using the 15, 25, 32 and 60 
Mq models using initial mass ranges of 12.5-20.0, 20.0-28.5, 
28.5-46 and 46-80 Mq in the IMF calculation. The overpro- 
duction factors for the CI - no20 case show a reduction in 
Sr isotopes to values just less than the 16 Ox2 line and a 
significant reduction in Y, Zr, Mo, Ru, Pd and Cd isotopes 
to values similar to the ST case and a reduction in Br and 
Rb isotopes to values close to the ls O/2 line. The branch- 
ing at 95 Zr is also affected, which mainly affects the relative 
overproduction factors of 96 Zr and 95 Mo. 



6 DISCUSSION 

The results in the previous section show that with an 
increased carbon burning rate, the contribution of the 
neutron-capture processes during hydrostatic burning stages 
to the yields of massive stars is modified significantly. 

The CU case exhibits a strong production of isotopes 
between the iron-group nuclides and the Ba-peak nuclides 



with regards to current massive star models (see Fig. 14 1 



This production originates from the s-process production in 
a convective carbon core in which mixing has caused the 
ashes of carbon burning to be transported out from the cen- 
tre of the star where it will be present in the supernova 
ejecta. This overlap was found in all but one of the CU 
models (15CU). Fig. 14 shows that the yields of the CU case 
are inconsistent with the weak s-process contribution to the 
solar system abundances (see for example the anomolously 
high abundance of Sr-Y-Zr peak and Ba-La peak nuclides 
compared to those with 60 < A < 90). Therefore, a strong 
resonance with (ury) — 6.8 x 10~ 5 eV at a centre-of-mass 
energy E com = 1.5 MeV in the 12 C + 12 C reaction rate is 
unlikely to be present in the reaction rate, according to the 
models used in the present analysis. 

For the CI case, an extended distribution is found but 
the overproduction factors are not as high as the CU case 



the Sr-Y-Zr peak and beyond, which is a signature dom- 
inated by the presence of overlap of a carbon shell with 
the convective carbon core. The large overproduction of Kr 
and Sr could suggest that the CI carbon burning rate is too 
high. In any case, it is unlikely that a solar metallicity model 
should demonstrate a strong overlap between the convective 
carbon core and the carbon-shell of the kind experienced in 
model 20CI. However, considering the present uncertainties 
in the stellar models such as the reaction rates (for example, 
the critical reactions 12 C(a,7) 16 and 22 Ne(a,n) 25 Mg), the 
initial composition and the treatment of convective-radiative 
boundaries, the abundance of Sr is not a significant enough 
constraint to assert that the CI rate would be inconsistent 
with the solar system abundance distribution. 

The production of Sr, Y, Zr and other heavier nuclides 
has been studied extensively as galactic chemical evolution 
models and observations have suggested the existence of an 
additional primary nucleosynthesis process, the lighter ele- 
ment primary process (LEPP) (Travaglio ct al. 2004; Montes| 



et al. 20071. The carbon core s-process and the mixing of 



heavy nuclei out from the centre could provide an alterna- 
tive nucleosynthesis scenario for the LEPP. It is tempting 
to underline the similarity between the LEPP signature and 
the anomalous carbon burning s-process component present 
in the CU models and partly in the CI models. However, 
we recall that the LEPP process should be primary if the 
solar LEPP and (low metallicity) stellar LEPP are indeed 
the same process (see for example Montes et al.|[2007 l. Al- 
though the carbon-core s-process features a primary neu- 
tron source, 13 C, the seed nuclei, 56 Fe, are secondary. Con- 
sequently, an s-process component using iron seeds in these 
conditions cannot reproduce the stellar LEPP abundances 
at low metallicity. Therefore, the carbon core s-process com- 
ponent is unlikely to represent the site for the stellar LEPP 
component at low metallicity. In addition, when the num- 
ber of seeds is lowered, the neutron captures per iron seed 



(see Fig. 15 I. The main nucleosynthesis differences occur at 



increases (see Eq. 4.1 1 and the distribution of s-process nu- 
clides extends to higher atomic mass. However, if the solar 
and stellar LEPPs differ in origin, the carbon core s-process 
may provide a solution to the solar LEPP. 



7 SUMMARY AND CONCLUSIONS 

In order to investigate the sensitivity of massive star evo- 
lution to the potentially large uncertainties in the carbon 
burning rate, fifteen stellar models with five initial masses 
of 15, 20, 25, 32 and 60Mq and three different carbon burn- 
ing rates were generated with GENEC and post-processed 
with the parallel post-processing code MPPNP. The yields 
for each model were then calculated and the consequences 
of the different rates on stellar evolution and nucleosynthe- 
sis were examined. The main conclusions are summarized as 
follows. 

An enhanced carbon burning rate directly affects the 
ignition conditions for carbon burning, which move to lower 
temperatures and densities. The reduced temperature low- 
ers the neutrino losses, causing the carbon burning stage to 
occur for a longer lifetime. An increasing dominance of neu- 
trinos formed through photoneutrino interactions is seen, 
rather than formation by pair-production. The change in 
temperature and the neutrino losses affect the convection 
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Figure 14. The overproduction factors of the predicted weak component for each rate, focussing on isotopes with atomic mass 50 < 
A < 150. Isotopes of the same element are connected by adjoining lines. The solid black line indicates the overproduction factor le O and 
the two dashed lines corresponds to the overproduction factors of ls O multiplied and divided by two. Changes to the overproduction 
factor of 16 O are negligibly small with changes to the carbon burning rate. The isotopic chains for Ge, Kr, Sr, Zr, Mo, Ru, Pd, Cd and 
Ba in the CU model are labelled for brevity. 
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Figure 15. The overproduction factors of the predicted weak component for the CI and CU rates relative to the ST rate. Isotopes are 
connected by adjoining lines. The isotopic chains for Kr, Sr, Zr, Mo, Ru, Pd, Cd and Ba in the CU model are labelled for brevity. 
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Figure 16. The overproduction factors of the predicted weak component relative to the solar system abundances with the 20Mq models 
removed from the calculation (CI - no20). The weak components for the ST and CI case including the 20 Mq models are included for 
comparison. Isotopes of a given element are connected by adjoining lines. The isotopic chains for Cu, Zn, Ge, Se, Kr, Sr, Y, Zr, Mo, Ru, 
Pd, Cd and Ba in the CU model are labelled for brevity. 



zone structure. In the models using the CI rate, the maxi- 
mum initial mass for the formation of a convective carbon 
core increases by a few solar masses from its current value of 
« 22Mq. In models using the CU rate, carbon core burning 
occurs in a convective core in the entire mass range. The in- 
creased carbon burning rates generally reduce the number of 
carbon burning shells (because they have a larger mass ex- 
tent) and increase the probability of overlap between differ- 
ent convective zones. Although the increased carbon burning 
rates used in this study strongly affect carbon burning, the 
impact on further burning stages (neon, oxygen and silicon) 
is small and does not present any clear trend. Therefore, 
no constraint can be applied to the 12 C + 12 C rate directly 
from stellar evolution considerations. 

The presence of a significant overlap between the con- 
vective carbon core and the convective carbon shell, as seen 
in most of the CU models and in model 20CI may present 
a further nucleosynthesis site worthy of investigation. This 
is especially true considering the present uncertainties in 
stellar models with regards to convective-radiative bound- 
aries and the abundance distribution exhibited by the car- 
bon core s-process. In particular, the carbon core s-process 
bares similarities to the solar LEPP. However, because of the 
secondary nature of the iron seeds, it cannot provide a solu- 
tion to the stellar LEPP at low metallicity. Further studies 
into the uncertainties relevant for low metallicity massive 
stars are required to confirm this statement. 

According to the present models, a strongly enhanced 
rate (the CU rate) due to the presence of a low energy res- 
onance (near to the Gamow peak) causes a large convective 



carbon core to exist in every stellar model. The large con- 
vective core will mix isotopes a considerable distance away 
from the centre of the star causing the ejecta to be polluted 
with matter rich in s-process isotopes. The overabundance 
distribution obtained with the CU rate is too high and has 
a vastly different shape. The yields are therefore incompati- 
ble with the weak s-process contribution to the solar system 
and the CU rate is therefore ruled out. 

A moderately enhanced rate (the CI rate), like the 
strongly enhanced rate, also affects the interior convection 
zones and consequently the structure of the star. With the 
CI rate, an overlap is only present in the 20Mq case, which 
enriches the ejecta with products of the carbon core s pro- 
cess. This enrichment predominantly involves nuclides at the 
Sr-Y-Zr peak and the heavier elements Mo, Ru, Pd and Cd. 
With this additional nucleosynthesis component, the over- 
production factor for Kr and Sr seems to be too high to be 
consistent with the solar system abundances since it would 
imply that the majority, if not all, of the solar Kr and Sr 
comes from massive stars, with only a smaller contribution 
from AGB stars at the Sr peak. For all the other masses, 
the changes in nucleosynthesis occur only from changes to 
carbon-shell burning, which are more subtle and involve iso- 
topes primarily at branching points. If the contribution from 
the 20 Mq model is not included (CI-no20), the yields ob- 
tained are very similar to the standard yields. Consequently, 
the CI rate is probably very close to the 'upper limit' for the 
carbon burning rate to lead to a weak s-process production 
compatible with the solar system composition. 

Given that an overlap between the convective carbon 
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core and shells has such a strong impact on the yields 
and that ID stellar models use the mixing length theory, 
which might not exactly represent the complex 3D nature 
of convective-radiative interfaces, it will be crucial to study 
such potential shell overlaps as well as overlap between burn- 



ing shells of different burning stages ( Arnett & Meakin|2011 



in 3D hydrodynamic simulations. It should also be acknowl- 
edged that the present conclusions are built on the assump- 
tion that the ratio of the a— and p-exit channels of the 12 C 
+ 12 C (13:7) reaction is preserved to lower energies. Further 
studies of this uncertainty, including also an analysis of the 
p-process in massive stars, will be discussed in a forthcoming 
paper (Pignatari et al. 2011, in prep.). 

The effects of the carbon burning rate on the stellar 
evolution and nucleosynthesis of massive stars demonstrates 
that nuclear physics experiments investigating the 12 C + 
12 C continue to remain relevant for the understanding of 
stars and further nuclear physics experiments, particularly 
at energies close to the Gamow peak for hydrostatic car- 
bon fusion, are highly desirable in order to improve stellar 
models. 
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APPENDIX A: PARALLEL-PROGRAMMING 
IMPLEMENTATION 

At a particular timestep, the parameters for a ID spheri- 
cal shell (or zone) are loaded into memory and a nuclear 
reaction network is calculated for that zone. This requires 
the inverse of a square matrix to be calculated, which has 
dimensions equal to the number of isotopes included in the 
network. For each timestep there are typically 10 3 zones, de- 
pendent on the stellar model and the evolutionary stage of 
the model, and there are ~ 10 6 timesteps per model. There- 
fore, the post-processing of a single stellar model requires 
~ 10 9 nuclear network calculations. With the nuclear re- 



action network specified in Table [T] including ~ 1.3 x 10 4 
reactions, the computational expense involved becomes sig- 
nificant; the typical duration of a single MPPNP run on a 
uniprocessor is approximately 10—12 months with current 
serial technology. Therefore, the application of parallel pro- 
gramming is an absolute necessity to allow the calculations 
to complete over a reasonable timescale. 

The choice of parallelism is a simple master-slave (or 
Workqueue) strategy where a single, master, processor al- 
locates work to a number of slave processors, which is im- 
plemented using the Message Passing Interface (MPI) li- 
brary routines in Fortran (Gropp et al. 19991. This is an 



implementation of parallelism where processors communi- 
cate information by passing 'messages' to each other with 
each processor having access to a local, private memory. The 
advantage of message passing is the ability to operate on 
distributed memory resources (such as cluster networks), as 
well as shared memory resources, and the ability to control 
explicitly how communications are handled and the parallel 
behaviour of the program. It is an embarrassingly parallel 
program^] which allows for an efficient parallelisation and re- 
duces dramatically the potential communication overhead. 
This was achieved by distributing 'work' over mass zones 
for each timestep, which are calculated independently from 
each other during the post-processing calculations. Here, a 
single unit of 'work' is defined as the nuclear reaction net- 
work calculation (in flops) for all involved species for a single 
zone at a particular timestep. 

The operation of the parallel program is as follows. 
First, the nuclear reaction rates and other global parameters 
are broadcasted to each slave so that each processor has the 
required data available in local memory. Then a loop over 
timesteps is entered. For each iteration of the loop a simple 
first-in first-out (FIFO) scheduler is invoked, which assigns 
work (in the form of a message containing the temperature, 
density and abundances) zone-by-zone (from the centre to 
the surface), first to all assigned processors and then to idle 
processors as they become available for further work. 

Load balancing is important to reduce the impact of 
idle processors on the performance. In MPPNP a simple 
load balancing scheme is specified, where the zones are allo- 
cated in order from the centre to the surface. This choice is 
made in lieu with the typical distribution of work over the 
interior of the star at any particular timestep. The distribu- 
tion is set by the dynamic network implemented in MPPNP, 
which adds or removes isotopes from the network calcula- 
tion depending on the nucleosynthesis flux limits (negligi- 
ble changes in abundances are ignored to save on unneces- 
sary computation). In general, the dynamic network assigns 
more isotopes to zones that have higher temperatures (since 
higher temperatures increase the nuclear reaction rates) and 
are convective (since the resultant mixing can cause an in- 
crease in the abundance of fuel). Therefore the distribution 
has a maximum in the centre and decreases with mass co- 
ordinate towards the surface, affected by the presence of 
convection zones. However, this is a general case; it is not 



6 An embarrassingly parallel program is one where slave proces- 
sors are not required to communicate information to each other 
during the run; the problem can simply be split and allocated in 
parts to a large number of processors. 
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unusual to have a non-monotonic distribution of work at 
particular steps throughout the model evolution, especially 
at the boundaries of convection zones and where neutron 
sources are efficient. 

The parallel burning step is followed by a serial mixing 
step. The change in mass fraction abundance of species i, 
Xi, over time, t, is calculated using the diffusion equation 



dX r 
dt 



d 
dm r 



no 2 n2 dXi 

D(47rr p) d^ r 



(Al) 



where m r is the mass coordinate (at radius r), p is the 
density and D is the diffusion coefficient calculated from 
mixing-length theory (MLT). The diffusion coefficient is nor- 
mally large enough (~ fO 16 cm 2 s _1 for hydrogen and he- 
lium burning) so that all convection zones, over a timestep 
At, act to smooth out immediately any sharp changes in 
abundance associated with concentrated nuclear burning. 

Figure |A1| shows the speed-up factor of MPPNP for a 
small test run (with 250 zones and 2000 timesteps; a typical 
stellar model uses ~ 10 3 zones and ~ 10 6 timesteps) com- 
pared to the theoretical laws predicted by Amdahl's and 
Gustafson's law for a program with a serial fraction of 1 per 
cent. Amdahl's law, 



S(P) = ^ = 



tp ft. + (1 - f)t./p 1+p-lf' 



(A2) 



gives the maximum speed-up, S(p), possible for a pro- 
gram with a fixed amount of work, i.e. the time spent run- 
ning serial computations is constant. In Eq. |A2| t s is the du- 
ration of the program with a serial fraction, /, on a unipro- 
cessor and t p is the parallel duration on a system with p pro- 
cessors. The close fit of this law with MPPNP suggests that 
the parallelisation is close to the ideal case and is not ham- 
pered by communication overhead or excessive initialisation. 
However, it would be preferable to achieve a parallelisation 
comparable to Gustafson's law, 



S(P) 



ft p +p(l-f)t p 



= P+/(1-P) (A3) 



which is the maximum speed-up possible with a con- 
straint on the parallel time, i.e. the time spent running par- 
allel computations is constant. This could be achieved by 
including more zones (for example, with the adaptive mesh 
refinement routine) , but the improved scaling would come at 
the expense of an increased workload. In any case, only 250 
zones were used in the test CcLSG j ELS the number of slave pro- 
cessors approaches 250, the total number of jobs allocated 
to each processor approaches unity. In this regime, the time 
spent by idle processors is likely to increase significantly and 
the speed-up factor will plateau. The post-processing calcu- 
lations for each model, using 60 slave processors, took ap- 
proximately 5 to 10 days each, depending on the model. 
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